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Abstract 

An elementary family of local Hamiltonians Ho/, I = 1,2,3,..., 
is described for a 2— dimensional quantum mechanical system of spin 
= I particles. On the torus, the ground state space Go/ is (log) 
extensively degenerate but should collapse under "perturbation" to 
an anyonic system with a complete mathematical description: the 
quantum double of the 50(3)— Chern-Simons modular functor at q = 
g27ri/£+2 -^yj^igj-^ -^g (,g^}j DEI. The Hamiltonian Ho/ defines a quantum 
loop gas. We argue that for ^ = 1 and 2, Go/ is unstable and the 
collapse to G^:/ = DEI can occur truly by perturbation. For £ > 3, 
Go/ is stable and in this case finding G^/ = DEi must require either 
e > > 0, help from finite system size, surface roughening (see 
section 3), or some other trick, hence the initial use of quotes " ".A 
hypothetical phase diagram is included in the introduction. 

The effect of perturbation is studied algebraically: the ground state 
space Go/ of Ho/ is described as a surface algebra and our ansatz is 
that perturbation should respect this structure yielding a perturbed 
ground state G^/ described by a quotient algebra. By classification, 
this implies G^/ — DEi. The fundamental point is that nonlinear 
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structures may be present on degenerate eigenspaces of an initial Ho 
which constrain the possible effective action of a perturbation. 

There is no reason to expect that a physical implementation of 
= DE^ as an anyonic system would require the low tempera- 
tures and time asymmetry intrinsic to Fractional Quantum Hall Ef- 
fect (FQHE) systems or rotating Bose-Einstein condensates — the 
currently known physical systems modelled by topological modular 
functors. A solid state realization of DE3, perhaps even one at a 
room temperature, might be found by building and studying systems, 
"quantum loop gases", whose main term is -ffo,3- This is a challenge 
for solid state physicists of the present decade. For £ > 3,i ^ 2 
mod 4, a physical implementation of DE£ would yield an inherently 
fault-tolerant universal quantum computer. But a warning must be 
posted, the theory at £ = 2 is not computationally universal and the 
first universal theory at ^ = 3 seems somewhat harder to locate be- 
cause of the stability of the corresponding loop gas. Does nature abhor 
a quantum computer? 
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Introduction 

In section 1, we write down a positive semidefinite local Hamiltonian i^o/ 
for a system of locally interacting Ising spins on a 2— dimensional triangular 
lattice or surface triangulation, i = 1, 2, 3, . . .. In the presence of topology, 
e.g. on a torus, Ho/ has a highly degenerate space Go/ of zero modes. On 
any closed surface Y, different from the 2— sphere, the degeneracy is polylog 
(2") = poly(w), where v is the number of sites in the triangulation and the 
2" is the dimension of the Hilbert space h of spins. On the torus the 
polynomial has degree = 1, when Y has genus g > 1 the polynomial has 
degree = 35^ — 3 (see Proposition 3.8). 

We argue for an ansatz (3.4) which exploits the peculiarly rigid algebraic 
structure of Go/ — it is a monoidal tensor category with a unique nontrivial 
ideal. The ansatz allows us to model any "perturbed" ground state space 
Ge,£ (which is itself stable to perturbation) uniquely as a known anyonic 
system or in mathematical parlance a modular functor . The functor is the 
Drinfeld double of the even — label — sector of the 5'[/(2)-Chern-Simons 
unitary topological modular functor at level £, DEi. Even labels corresponds 
in physical terms to the integer spin representations so the even-label-sector 
derives from the group 5*0(3). 

The Hamiltonian Ho/ defines a quantum loop gas which can be compared 
(see Figure 3.12) with the classical analog. The statistical mechanics of 
classical loop gases [Ni] identifies a known critical regime and from this we 
infer that for i = 1 and 2, Go/ is unstable and the collapse to G^/ = DEi 
is truly by perturbation, for i > 3, Go/ is stable and in this case finding 
G^/ = DE£ requires e > > 0, or some other device (see section 3), hence 
the initial use of quotes " ". Figure 0.1 is a hypothetical phase diagram. 
The stability of Go/ at £ = 3 is probably very slight — see footnote 6 in 
section 3 and the corresponding discussion. 
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The reader should not be alarmed that a "doubled" Chern-Simons theory 
arises. The doubled structure makes it a gauge theory and, as we will explain, 
the double, being achiral, is more likely to have a robust physical realization. 
The modular functor DEi has A = ( J ) "labels" or , physically, A super 
selection sectors for quasiparticle excitations (including the empty particle.) 
Physically this means that a local bit of material, a two dimensional disk 
with a fixed boundary condition, which is in its unique ground state G^^^ 
can have A types of point-like anyonic excitations (presumably with expo- 
nentially decaying tails) which can only be created in pairs. A is the number 
of order integers pair {x, y) with < x,y < i, and x,y — even. By math- 
ematically deleting small neighborhoods of such excitations a ground state 
vector is approximately achieved in the highly degenerate ground state space 
associated to a punctured disk with boundary conditions. It is known 
[FLW2 ] and [FKLW] that for £ > 3, £ 7^ mod 4, there is a universal, in- 
herently fault-tolerant, model for quantum computation based on the ability 
to create, braid, fuse, and finally distinguish these excitations types. Thus 
Ho £ could be of technological importance: a physical system, a "qiiantum 
loop gas" , in this (perturbed) universality class could be the substrate of a 
universal fault tolerant quantum computer. 

Any unit vector ^ e Go,e is a superposition of classical ±— spin states |^) 
which are distinguished by the eigenvalues ±1 of a commuting family of Pauli 
matrices equal to | J -1 1 at vertex v. Sampling = Sai|\t'j) by observing 
{crl}, we "observe" a classical with probability |ajp. The domain wall 7j 
separating the -|- - spin regions from — spin regions of may be thought of 
as a random, self dual, loop gas [Ni]. This random state is self dual because 
there is a symmetry between "up" and "down". It is a Gibbs state with 
parameters k — 0, self dual, and n = (2 cos ^) , where the weight of a 
configuration 7 is u;(7) = e-*^^*"*^^ ^^"^th 7)^* components 7, it jg known that for 
< n < 2 and k = the loop gas is critical, sitting at a 2"^^ order phase 
transition as k crosses from negative to positive. This information together 
with sections 3 or 4 support a phase diagram like the one shown in Figure 
0.1 with parameters d :— 2cos;^ = y/n and e. The parameter e scales a 
local perturbation term eV. We will argue that the simplest choice for V, 

V — ( E (7* ) , is a likely candidate. The diagram is labelled "hypothetical" 

\site i / 

since there is no proof of its accuracy. 

The challenge to solid state physics is to find or engineer a two dimensional 
quantum medium in the universality class, DE3 below. 
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hypothetical phase diagram for ground state space 
G,, of < id] 

\ sites V ) 

Figure 0.1: Shaded regions are the topological phases 

DEI, DE2, DE3, DEA, . . .. Doubly shaded regions are the computationally 
universal phases DE3, DE5, DE7, . . .. We have no way of predicting if the 
topological phases are actually in contact with each other as drawn. Solid 
lines are phase boundaries between inequivalent systems. 
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Figure 0.2: For sufficiently many particles we have recorded if the (general- 
ized) braid group representations are: T (trivial), A(abelian), N(nonabelian), 
U (computationally universal), we have called the total number of elementary 
particles, including those that reverse the (|+), |— )) coloring, "specific heat" 
as it counts local degrees of freedom above the ground state. The color con- 
stant elementary particles arc the irreducible representations of the corre- 
sponding linear category (see § 2). Coloring-reversing particles are explained 
at the end of ^ 2. 
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The presumptive approach — nearly universal in the literature — to build- 
ing a quantum computer is quite different from our topological/anyonic start- 
ing point. It is based on manipulating and protecting strictly local — as 
opposed to global or topological — degrees of freedom. It may be called the 
"qubit approach" since often a union of 2— level systems (with state space 
(^Cf ) is proposed. Actually, the number of levels — or even their finiteness 

i 

— is not the essential feature, it is that each tensor factor of the state space 

— call it a qunit — is physically localized in space (or momentum space). 
The environment will - despite the best efforts of the experimentalist — in- 
teracted directly with these "raw" qunits. It has long been recognized ([S], 
[U]) that the raw qunits must encrypt fewer "logical qubits." The demon in 
this approach is that very low initial (or raw) error rate - perhaps one error 
per 10~^ operations - and large ratios of raw to logical qunits ~ 10^ seem 
to be required [Pr], to have a stable computational scheme. This problem 
pervades all approaches based on local or "qubit" systems: liquid NMR, solid 
state NMR, electron spin, quantum dot, optical cavity, ion trap, etc. . .. 

Kitaev's seminal paper [Kl] on anyonic computation, amplified in numer- 
ous private conversations, provides the foundation for the approach described 
here. Anyons are a (2 + 1)— dimensional phenomena: when sites containing 
identical particles in a 2— dimensional system are exchanged (without colli- 
sion) there are, up to deformation, two basic exchanges; a clockwise and a 
counter-clockwise half turn - or "braid" if the motion is considered as gen- 
erating world lines in 2 -|- 1— dimensional space-time. The two are inverse 
to each other but of infinite order rather than order = 2. So whereas only 
the permutation needs to be recorded for exchanges in E?, in "statistics" 
becomes a representation p of a braid group B into the unitary group of 
a Hilbert space h encoding the internal degrees of freedom of the particle 
system: 

p:B — > U{h). 

Since a representation into unitary transformations, "gate set" in (quantum) 
computer science language, is the heart of quantum computation it is not 
really a surprise that any kind of particle system with a sufficiently general 
(it certainly must be nonabelian) image p{B) can be used to build a universal 
model for quantum computation. This has been shown in [FLWl], [FLW2], 
and [FKLW]. 

What are the advantages and disadvantages of anyonic verses qubit com- 
putation? The most glaring disadvantage of anyons is that no one is abso- 
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lutely sure that nonabelian anyons exist in any physical system. Two dimen- 
sional electron liquids exhibiting the fractional quantum Hall effect FQHE 
are the most widely studied candidates for anyonic systems. The Laughlin 
state at filling fraction v = 1/3 has observed excitations charges of (l/3)e 
and these are convincingly linked by the mathematical model with a sta- 
tistical factor of a; = ie^'^*/^ for the exchange of such pairs. Quasiparticle 
excitations with nonabelian statistics is one of the most exciting predictions 
of Chern-Simons theory as a model for the FQHE. With a few low level (e.g. 
£= 1,2 or 4, when G = SU{2)) exceptions nonabelian anyonic systems are 
capable, under braiding, of realizing universal quantum computation [FLW2]. 
The essential point is that the "Jones representation" of the braid group (on 
sufficiently many strands) associated to the Lie group SU (2) has a dense 
image at least in SU{h) C U{h), h an irreducible summand of the represen- 
tation. At u = 5/2 according to [RR] the Hall fiuid is modelled by a U{1) 
theory coupled to CS2 [the Chern-Simons theory of SU{2) at the 4**^ root 
of unity (level i — 2)]; the latter is a theory with a nonabelian "Clifford 
group" representation. This model was selected from conformal field theo- 
ries to match expected ground state degeneracies and central charge, and is 
further supported by numerical evidence on the overlap of trial wave func- 
tions. Though very interesting, this representation is still discrete and is not 
universal in the sense of [FLWl]. However at = 8/5, with perhaps weaker 
numerical support [RR], it is thought that the Hall fiuid model contains by 
CSS (level= 3, 5*^^ root of unity). Here braiding and fusing the excitation 
would yield universal quantum computation [FLWl]. 

So let us, for the sake of discussion, accept that FQHE systems have 
computationally universal anyons, we are still a long way from building a 
quantum computer. FQHE systems are very delicate: 

1. The required crystals have been grown successfully only in a few labo- 
ratories. 

2. The temperatures at which the finer plateaus are stable are order milK. 

3. The chiral asymmetry intrinsic (For CS2 and CS3 the central charge 
is I and | respectively.) to the effect requires an enormous transverse 
magnetic field, order 10 - 15 Tesla to reduce magnetic length to where 
conduction plateaus are observed. 
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At feasible magnet lengths]]], the Coulomb interaction between electrons is 
at least three orders of magnitude weaker than in solids. Corresponding to 
the weakness of these interactions the spectral gap protecting topological 
phases is necessarily quite small. Perhaps for this reason, even the most 
basic experiments to prove existence of "nonabelions" have not been carried 
out, and the use of these systems for computations appears unrealistic. 

For applications such as breaking the cryptographic scheme RSA, it can 
be estimated that several thousand anyons must be formed, braided at will 
(perhaps implementing tens of thousands of half twists), and finally fused. 
This appears to be a nearly impossible task in a FQHE system. 

The main point of this paper is that computationally universal anyons 
may be available in more convenient systems. Ho/ is a local model for a 
paramagnetic system of Ising spins with short range antiferromagnetic prop- 
erties. Written out in products of Pauli matrices Ho/ is seventh order (on 
the standard triangular lattice) and thus looks complicated compared to, say, 
the Heisenberg magnet. But geometrically it is quite simple and its ground 
states are known exactly. A 2— dimensional material in the universality class 
DEI proposed as the ground state space for H^ i,i ^ 1,2, or = mod 4, 
will have excitations - "quasiparticles" - capable of universal fault tolerant 
quantum computation within a model that allows creation, braiding, fusion 
and measurement of quasiparticle type. 

A topological feature is not too easy to detect; by definition, topological 
properties cannot be altered or measured by purely local operators but in- 
stead require something akin to an Aharonov-Bohm holonomy experiment. 
So perhaps the universality class of H^/ already exists in surface layer physics 
but is waiting to be discovered. Or perhaps with Ho/ in mind something in 
its (perturbed) universality class can be engineered. If this is possible there 
would be no reason to expect the system to be particularly delicate. The 
characteristic energies for magnets are often several hundred Kelvin [NS]. 
Furthermore the modular functor DEi (this includes the information of the 
various braid group representations, 6j— symbols, S and fusion matrices) 
which arises is amphichiral, the central charge c = 0, so there is no rea- 
son that time symmetry must be broken and no apparent need for a strong 
transverse magnetic field. These two features are in marked contrast to the 

-'^In semiconductors with dialectic constant e « 10 and \B\ k, 10 Tesla the characteristic 
length I = yJeh/eP ~ 150A compared to about AA separation between the ions in a 
crystalline solid. 
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delicate FQHE systems. 

Subsequent to the initial draft of this paper a different local Hamilto- 
nian H'^ ^ was found which bears the same relation as ifo/ to the topological 
modular functors DEi, but has potential advantages: 

1. it is expressible as 4*"^ rather than 7^^ order interactions and 

2. its classical analog is the much studied self dual Potts model for q = 
(2cos^)^ 

We have added a section 1' following section 1 to explain this alternative 
microscopic model. 

We make no proposal here for a specific implementation of Ho/ or for how 
to trap and braid its excitations but we hope that models in the spirit of [NS] 
for the high Tc cuprates may soon be proposed. In this regard, we note that 
relatively simple - but still non classical-braiding statics have been proposed 
[SF] in conjunction with the phenomena of spin-charge separation [A] for high 
Tc cuprate super conductors above their Tc. Certain —1 phases are predicted 
to occur when braiding the electron fragments "visions" and "chargeons" 
around each other and around ground state defects called "holons". Also 
contained in this paper is the suggestion that topological charges might in 
passing though a phase transition become classical observables, e.g. mag- 
netic vortices. Similarly other phase transitions might link higher {£ = 3) 
to lower {£ = 1) topological phases and might be useful in measuring quasi- 
particle types. Whether even the simplest topological theory is realized in 
any known superconductor is open, but [SF] is cited as precedent for anyonic 
models for solid state magnetic systems with high characteristic energies. So 
while the FQHE motivates this paper, we hope we have steered toward its 
mathematical beauty and away from it experimental difficulties. 

What are the generic advantages of anyonic computation? First infor- 
mation is stored in topological properties "large scale entanglement" of the 
system that cannot be altered (or read) by local interaction. This affords a 
kind of physical stability against error rather than the kind of combinatorial 
error correction scheme envisioned in the qubit models - "hardware" rather 
than "software" error correction. Second, at least in the simplest analysi^, 

^Kivelson and Sandih[KS] find that Landau level-mixing in FQHE can thicken the tails 
to polynumerical decay, but this is not a fundamental effect. 
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one expects excitation of a stable system to be well localized with exponen- 
tially decaying tails. Thus physical braiding should approximate mathemat- 
ical braiding, p : B — > U{h) up to a "tunnelling" error of the form e~'^^, 
where c is a positive constant, and L is a microscopic length scale describing 
how well separated the excitations are kept during the braiding process. This 
error scaling is highly desirable and seems to have no analog in qubit mod- 
els. While tunnelling treats virtual errors, errors which borrow energy briefly 
from the vacuum, actual errors would be expected scale like e^^'^'^/'^° where 
To is a character energy for the system and T the operating temperature. 
This is essentially the error analysis Kitaev made for his anyonic system, the 
toric code [Kl]. 

This paper draws on three sources of inspiration: 1) Kitaev's paper [Kl] 
on anyonic computation, 2) the FQHE, and 3) rigidity in the classification of 
von Neumann algebras subfactors. Rigidity implies that certain monoidal 
tensor categories have very few ideals. But when interpreted physically, 
"ideal" means "definable by local conditions" , so we find that a certain lo- 
cality assumption (Ansatz 3.4) strongly limits the physics. This provides an 
algebraic approach to the perturbation theory of Ho/ - and perhaps yields 
greater insight than would be possible by analytic methods. We find that 
for Ho/ the polylog extensively degenerate space of 0— modes Go,e possess 
in addition to its linear structure, an important "multiplicative" structure 
— the structure of a monoidal tensor category - which we argue, should be 
preserved under a perturbation. The rigidity of type IIi factor pairs, an as- 
pect of which is stated as Thm 2.1, provides a unique candidate for the (still 
finitely degenerate) "perturbed" ground state space G^/ of if^,^. The space 
is a braid group representation space with the representation induced by 
an adiabatic motion of quasiparticle excitations. 

Throughout, the excitations on a surface Y are assumed to be localized 
near points so excited states of H^/{Y) become ground states of H^/{Y~) but 
now on a punctured surface Y~ with "boundary conditions" or more exactly 
"labels," (see section 2.) We treat excited states indirectly as ground states 
on the more complicated surface Y~ . 

The existence of a stable phase G^/ = DEI will be argued by anal- 
ogy with the FQHE where topological phases are found to be stable, from 
algebraic uniqueness, and via "consistency checks". But these arguments 
constitute neither a mathematical proof nor a numerical verification. The 
latter may be exactly as far off as a working quantum computer. It was 
precisely the problem of studying quantum mechanical Hamiltonians in the 
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thermodynamic limit, e.g. questions of spectral gap, that lead Feynmann 
[Fe] to dream of the quantum computer in the first place. It is curiously self- 
referential that we may need a quantum computer to "prove" numerically 
that a given physical system works like one. 

We turn now to the definition of Ho,e and H'^ f, returning later to am- 
plify on the relations to quantum computing, C*— algebras, Chern-Simons 
theory, and topology. (The connection between Chern-Simons Theory and 
complexity classes is discussed in [Fl].) 

In addition to Alexei Kitaev, I would like to thank Christian Borgs, Jen- 
nifer Chayes, Chetan Nayak, Oded Schramm, Kevin Walker, and Zhenghan 
Wang for stimulating conversations on the proposed model. 

1 The model 

The model describes a system of spin = | particles located at the vertices v 

n 

of a triangulated surface Y. The Hilbert space is 7:^ = where is the 

v=l 

local degree of freedom {| -|-), | — )} at the vertex v. The basic Hamiltonian 
Ho/ is written out below as a sum of local projections and thus is positive 
semidefinite. The ground state space (energy = vectors) Go/ of Ho/ can 
be completely understood (this is unusual since the projector above do not 
commute) and identified (as n — > oo) with what we call the even Temperley- 
Lieb surface "algebra" ETL^ where d — 2cos;^. 

Ultimately our focus will be on the ground states on a multiply punctiued 
disk — the puncture corresponding to anyonic excitations (see section 5). 
Two issues arise: (1) non-trivial topology and (2) boundary conditions. The 
boundary conditions are quite tricky so it is best to work first with closed 
surfaces of arbitrary genus (even though these are not our chief interest) 
to understand the infiuence to topology alone "liberated" from boundary 
conditions. 

Y will denote a compact oriented surface throughout. In combinatorial 
contexts, Y will be given a triangulation A with dual cellulation C. Initially, 
we consider the case where Y is closed, boundary Y — dY — 0. We will speak 
in terms of the dual cellulation by 2— cells or "plaques" c. For example, if 
y is a torus it may be cellulated with regular hexagons. This is a perfectly 
good example to keep in mind but higher genus surfaces are also interesting, 
while the sphere is less so. Soon we will consider surfaces with boundary. 
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Distributing over 0, one writes Ti = span {classical spin configu- 
rations on plaques} =: span {sj}. Let c be a plaquet, s a classical spin 
configuration and s*^ that configuration with reversed spin (+ — > — and 

> +) at c. For 1 < i,j < 2" define hij{c) = 1 if {l)sj = s^ and (2)sj 

assigns the same spin ± to c and all its immediate neighbors, and hij{c) — 
otherwise. Define gij{c) — 1 ii {l)sj — and (2) the domain wall 75. between 
+ and — plaques, in the spin configuration Sj, meets c in a single connected 
topological arc, and gij{c) = otherwise. We define: 



Ho,e = Yl 3ij{c) {{\si) - \sj)) {{si\ - {sj\)) + 

plaques c, pairs 
of spin states Si,Sj 

plaques c, pairs ^ ^ ' ^ 

of spin states Si,Sj 

The constant k is positive and may, in this paper, be set as k = 1. To 
help digest the notation each of the two sums has n2^" terms most of which 
are zero. It is easy to see that gij — gji. If the domain wall 7 meets c in 
a topological arc reversing the spin of c isotopes the domain wall across c 
to the complementary arc = 5c \ 7. Contrariwise if hij = 1 then hji — 0. 
The parameter d could be any positive real number but we will be concerned 
mainly with d = 2 cos = 1,2,3, The cases i = 2, d = \/2 and 

£ = 3, d = ^-^y^, the golden ratio, are of particular interest. Finally, each 
term in the definition of Ho/ should be read, according to the usual ket- 
bra notation, as orthogonal projection onto the indicated vector: — \sj) 
or |sj) — l/d\sj). These vectors (whose projectors occur nontrivially in the 
sums) are certainly not orthogonal to each other (using the inner product |+) 
hermitian orthonormal to |— ) in C^, extended to define the tensor product 
Hermitian structure on H) so those individual projectors do not commute. 
It is therefore surprising at first that we can completely describe the (space 
of) zero modes Go/ of this positive semidefinite form. Ho/. However once 
the description is given the surprise will evaporate for it will be clear how 
Ho/ was "engineered" precisely to yield this result. Identifying Go/ is the 
next goal. 

Associate to the closed oriented surface Y an infinite dimensional vector 
space ETLd(y), the even Temperley-Lieb space of Y. It is the C— span 
of "isotopy classes" of closed bounding 1— manifolds 7 modulo a relation 
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called d— isotopy. The "bounding" condition means that 7 is a domain wall 
separating Y into two regions which could be labelled "|+)" and a "|— )". 
Neither 7 nor the regions are presumed to be connected. We do not orient 
7, so we do not distinguish here between states which differ by globally 
interchanging |+) and |— ). The term "1— manifold" means 7 does not branch 
or terminate at any point. Isotopy, of course, means gradual deformation. 
The (i— isotopy relation: 7 — (i(7\7o), when imposed, says that if a component 
7o of 7 bonds a disk in Y then 7 = ^(7 \ 70), d times the value on the 
submanifold with 70 deleted. We often work with the dual ETL*^{Y), which 
are the functions / on bounding isotopy classes satisfying 7(7) = (i(/(7\7o))- 
Let 7' be a 7 as above enhanced by one of the two choices for "signing" 
the complementary regions. Define ELT\{Y) to be C— span {7'}, so that 
ELT^^iY) are the functions from {7'} obeying the rf— isotopy relation. 

Both the definition of i^o/ and ETL*^ can easily be extended to F a 
compact surface with boundary = dY, given a fixed boundary condition, 
the points where 7 meets dY (transversely). So if A is a triangulation of 
(y, dY) with dual cellulation C and if the spin configuration |+) or | — ) is fixed 
at every vertex (= dual cell) on dY then formula (1.1) defines a Hamilto- 
nian operator on the configurations with that boundary condition provided, 
in both terms, we restrict the sum to plaques c which do not meet dY. 
This prevents "fluctuations" from altering the boundary conditions. Deflne 
Ho/{Y, dY) in this way. Similarly if a 2— coloring (or +, — "signing") of Y is 
fixed along dY we may consider a relative 7' as a extension of this signing 
to a division of Y into + and — signed regions (which are presumed to lie 
in y as subsurfaces). Now relative to the boundary condition (the signing) 
ETV^(Y, dY) is defined as functions from {7'} to C which obey the d— iso- 
topy relation; ETL*^{Y, dY) is the set of such function invariant under ~, the 
global 1+) < > |— ) interchange. 

If A is a triangulation on a surface Y, with or without boundary then we 
have the combinatorial versions of ETL|^(F) and ETLrf(y), ETLJ'(F) and 
ETL^(y) (resp.) define using only (1+), |—)) 2— colorings in which each dual 
2— cell (plaquet) is + or — . 

There are natural maps of C— vector spaces: 

ETLf(y) ^ ETL)i(y) and ETLj(y) ^ ETLrf(y). (1.2) 

These maps are of course never onto (only the simpler rf— isotopy classes are 
realized). Also for certain triangulations A the kernel can also be non-zero 
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(due to "stuck" configurations). However, it is easy to see tliat as A is 
subdivided. ETL^(Y) approximates ETLaiY) in the sense tiiat the direct 
Umit \jmETL^{Y) = ETLd{Y), similarly \uaETL^'{Y) = ETL-^{Y). 

Let A be a fixed triangulation of Y (with fixed boundary condition — 
a ^ — projective (|+), |—)) 2— coloring — if dY ^ 0), set Go/(^, A) — zero 
modes (ground state space) of the positive semidefinate ifo,^ defined above 
(1.1). Clearly H^^ is "—invariant and so Go/ is "—invariant. Note that 
is not always fixed point free: on F = T^, the configuration which is |+) on 
an essential annulus A C and |— ) on \ ^4 is a ~ — fixed point. Let 
GjJ'^(y, A) denote the +1— eigenspace of ~. 

Proposition 1.1. For Y a closed surface or a surface with fixed boundary 
conditions, there are natural isomorphisms Go,e{Y, A) = ETL^'(y) and G^^{Y, A) = 

ETLj(r). 

Proof: From line (1.1), * G Go(F, A) iff ^ = E{f\si))\si) for some linear 

i 

functional obeying the d— isoptopy relation, thus Go/{Y,A) = ETL^'(y). 
The involution ~ acts compatibility on both sides so ETL^ (F) may be iden- 
tified as the +1 eigenspace of ~ on the r.h.s. □ 
When we come (section 3) to imposing the mathematical structure of 
a modular functor (or TQFT) on the ground state spaces G^^iY) for var- 
ious surfaces Y we will need to impose a base point * on each boundary 
component C C dY. This is directly analogous to the framing of Wil- 
son loop in [Wi], in fact the base point moving in time defines the first 
direction of a normal frame to the Wilson loop in the 2-1-1 dimensional 
space-time picture. As in the previous application, the base point is intro- 
duced for mathematical rather than physical reasons. It allows the state 
vectors in each conformal block to be identified precisely and not merely 
up to a (block-dependent) phase ambiguity. Concretely in our model the 
base point prevents domain walls from spinning around a puncture. Note 
that if (a superposition of) domain walls 7 represent an eigenspace for Dehn 
twist around the puncture with eigenvalue A 7^ 1 and if twisting is not pre- 
vented then the relation I7) = AI7) will occur, killing the state I7) which is 
certainly not desired. I thank Nayak for pointing out that although choos- 
ing base points breaks symmetry, none of the physics depends on which 
base points are chosen. The Hamiltonian has a [/(I) x • • • x C/(l)— gauge 

k 

symmetry where k — ^ boundary components of Y. 
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1' An alternative microscopic model. 

In this subsection we present an alternative Hamiltonian, H'^ ^, on a cellulated 
surface {¥,€). We do not restrict to triangulation since the square lattice 
actually yields the simplest form. It has the same relation, in the infrared, 
to topological theories as does Ho/. In this model the degrees of freedom 
are on bonds and the loops lie in a "midlattice" separating the |+) clusters 
from the |— ) dual clusters (isolated vertices and isolated dual vertices count 
as clusters). There are perhaps three advantages: 

1. On a square lattice, all terms in the Hamiltonian have order 4 (as com- 
pared to seven in the previous model). It is simple enough that we 
expand it ClS db product of Pauli matrices. 

2. The corresponding classical statistical mechanical model is the Potts 
model in cluster expansion (FK) form with q = (2cos^)^.; and 

3. The loops in this model are "fully packed" so no isotopy is possible, 
only (i— isotopy. In particular the total length of the loops separating 
1+) from |— ) is configuration independent. Here is H'^^^ ; the notation 
is explained below. 



K is a positive constant which for symmetry we suppose to be k = 1. 
Again d = 2cos^. On each bond there is a spin = ^ degree of free- 
dom = = span {|+), |— )}. The first summation is over all plaques 
(2— cells) with a market edge (So, if the surface is a 10 x 10 torus cellu- 
lated with 100 squares, the first summand contains 400 terms.) Each term 
in the first sum is orthogonal projection onto the vector 
(|-) ® 1+) ® 1+) (S) ■ ■ ■ ® 1+) - ^1+) ® 1+) ® 1+) ® ■ ■ ■ ® 1+)) where the ten- 
sor factors begin with the bond containing the dot and proceed counterclock- 
wise around the plaquet. (Of course this projector is understood to be ten- 
sored with the identity over all remaining bonds.) Perhaps this is confusing. 



^:.-E(|3)-^|4)) ((3|->i) + 



(1.3) 
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but we have used the notation |3) for the first and |4) for the second basis 
vector in this combination because in the square lattice case, those numbers 
count the + signs: a more elaborate notation would be \ni — 1) — 
The second term is the "double dual" of the first where one duality swaps 
cellulation with dual cellulation (homology with cohomology) and the other 
duality swaps |+) and |— ). Thus the second summation is over vertices with 
a marked incoming bond; the vector |1) denotes |+) (8) |— ) (8) |— ) ® • • • ® |— ) 
and |0) denotes |— ) C?) |— ) C?) |— ) ® ■ ■ ■ (E) |— ), again reading counterclockwise 
from the dot. The ^ reminds us that we are reading around a site not a 
plaquet. In the case of the square lattice the two types of terms may be 



1th 



expressed 

rr - I 1 I 
— I 1 I- 

£ has two types of terms: 



degree polynomial in Pauli matrices: cr^ 



l+> I-) 

1+) I 1 1 
l-> I -1 I 



and 



1 

16 

1 
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d'' ' d^ 



® [(/ + al) 0(1 + al) ® (7 + erf)] , and 
® [(/ - al) (/ - al) ® (/ - al)] 
[{I + al)0{I + al)0{I + al)], and 
[(/ - al) ® (/ - al) ® (/ - al)] 



The proper context for understanding H'^ ^ is Baxter's "mid lattice" [B] . 
If C, C* are cellulation and dual cellulation, let d — cH c* he the general 
intersection of a plaquet and a dual plaquet. Put a center * in c' and a center 
point in each of its boundary 1— cells. Subdivide d by the cone 

n 

of IJ *j to * and let c" denote the general plaquet of this subdivision. The 

i=l 

collection {c"} are precisely the plaquets of the "mid lattice" . As an example, 
for the square lattice of unit size the dual lattice is shifted by (1/2, 1/2) and 
the resulting mid lattice is spanned by vectors {(0, 1/4), (1/4,0)}. 

A classical configuration s : {bonds of C} — > \ ~)} is encoded as the 
union of bonds on which s is |+), the components of which are called clusters 
and the union of the duals of bonds on which s is |— ), whose components 
are called dual clusters. There is a well defined 1— manifold (multi— loop) 
7s in the mid lattice which separates clusters from dual clusters. 
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The Hamiltonian H'^^ builds in dynamics which fluctuates broken (|3)) 
and complete (|4)) boxes and broken |1) and complete (|0)) dual boxes with 
a prescribed weight factor = d. The vector |4) encodes a small face— centered 
loop, 0/ in the mid lattice while |0) encodes a small vertex— centered loop, 
Oy. The first term projectors, by annihilating |3) + (i|4), enforce a relation . 
If ^ = J2i O'i^i £ G'^,e: the zero modes of if^^, and i is written out as (index 
on boundary plaquet, distant indices k) then, 

da,„. — a... ^, for all k. (1.4) 

Examining this relation on mid lattice multiloops 7 (and suppressing k) we 
see that 7|4) differs from 713) in that an Of has been added to the isotopy 
class of 7|3) by "pinching off' a small bend in 7|3). Correspondingly 7|4) has 
it's coefficient a|4) equal to d times the coefficient a\3) of 7|3). Similarly for the 
double dual: up to isotopy 7|q^ = 7|j;^ UOj; and for a zero mode the coefficients 
much satisfy: 

da\i),j, = a\d),n (1-5) 
analogous to line (1.2) and Proposition 1.1 we have: 

Proposition 1.2. There are natural maps: ETL'^'iY) — > ETV^(Y) and 
ETL'^iY) — > ETLd{Y). In the (direct) limit they become isomophisms. 
There are natural isomorphisms: G'^^{Y,C) ^ ETL'^J{Y) andC'^^^ ^ ETL^^{Y) 

□ 

Proposition 1.2 replaces the triangulation A with the cellulation C, so 
ETLf^- means formal configurations, Ea^s, where s assigns |+) or |— ) to 
the bonds of C and obeys (1.4) and (1.5). ETL^ are formal configurations 

which are also invariant under the global swap, (~), |+) < > |— ). In the limit 

this relation expresses d— isotopy of the mid lattice domain wall 7. Note, 
however that the first two maps mentioned in the proposition 1.2 are not 
necessarily injective. The situation is summed up by the following example. 
On a 2 X 2 square torus the two possible staircase diagonals i.e. |+) on 
one positively sloping diagonal and |— ) on the complement, to not fluctuate 
(are not in the same ergodic component) whereas already in the 3x3 torus 
there is enough room that any two staircases of slopes = 1 are connected by 
fluctuations. 
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, however 



= d 



by symmeti'y 



Figure 1.3 



Remark 1.3. Because of their importance in solid state physics, we observe 
that a certain ring exchange Hamiltonian H'^ is the parent of all H'^ ^ in that 
the zero— modes G'^ contain the zero modes G'^^, for all I. Each G'^^ arises 
from a distinct linear constraint on G". 



= E - |3')) ((3| - (31) + - - 

□ ■+ 

1 3') is like |3) except cycled one step: 

|3') = 1+) (g) |-) ® 1+) ® . . . , similarly 

|T') = |-)®|+)® |-)®---®|-). 

Note that (|3) - |3')) G span ((|3) - i|4)) , |3') - i|4)), etc., so G;_^ C G'J. 
The zero modes G" can be identified with the (combinatorial) isotopy classes 
of domain walls between |+) and |— ) regions. 

Measuring spins (by a family of a^'s) converts a ground state vector \1/ G 
Go/ or G'o^ into a classical probabilistic state = meas.(\l/) which turns out 
to be a Gibbs state. The statistical physics of meas.(\l/) plays an important 
role in section 3. First, however, we use section 2 to lay down the algebraic 
framework. 



2 Things Temper ley- Lieb 

The generic Temperley-Lieb algebra is a tensor algebra over the complex 
numbers adjoined an indeterminate d. Often d is written in terms of an- 
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other indeterminate A as d = —A^ — A^'^ . In degree n the algebra TL„ has 
generators 1, ei, . . . , e„_i and the relations ef = Cj, CjCj = ejCj if |i — j| > 2 
and ejCjiiej = ^e,. Pictorially, after V. Jones and L. Kauffman, we may 
think of the generators as pictures of arcs disjointly imbedded in rectangles 
(multiplied by the coefficient l/d) and multiplication as vertical stacking. 
For example for n = 4, we have: 



U 

n 



u 
n 



u I 

n : 



Figure 2.1 



There is a convention that any closed circle ( and these may arise when the 
pictures are stacked) should be regarded as a factor of d. All closed circles in 
a picture should be deleted and the resulting picture should then be formally 
multiplied by d'^'^™'^^^^\ The reader can now easily verify the relation by 
stacking pictures. Kauffman proved the algebra of such pictures has no other 
relations [K] . A tensor structure between grades TL„ ^ TL^ — > TL„+m is 
created by horizontal stacking. An inclusion TL„ — > TLn+k is obtained by 
adding k vertical strands on the right. The union of grades is the (generic) 

oo 

Temperley-Lieb algebra, TL = U TL„. The structure of this algebra is 

ra=l 

completely worked out in [J]: Each grade TL„ has dim(TL„) = the 
v}"^ Catalan number, and is a direct sum of matrix algebras that fit together 
via a rather simple Brattelli diagram. Also of interest are specializations 
where the indeterminate d is set to a fixed nonzero real number. Here the 
structure differs from the generic case only when d assumes a "special value" 
d = 2 cos £ a positive integer, and has been worked out by Goodman and 
Wenzl [GW]. 

There is an involution - on TL which acts by refiecting the rectangle in a 
horizontal line and conjugating coefficients {d = d) making TL a * - algebra. 
Using this, the "Markov trace pairing < a,b >:= trace (ab) may be defined. 
"Trace," on pictures, means closing a rectangular diagram by a family of arcs 
sweeping from top to bottom and then evaluating each circle as a factor of d 
times 1 G C. Extend this definition to a Hermitian pairing on TL 
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Figure 2.2 



Theorem 2.1. ([J]) The trace pairing (, ) : TL(^TL — > C[(i], when d is 
specialized, to a positive real number, becomes a positive definite Hermitian 
pairing ( , )d ■ TLrf(^TL(i — > C exactly for d > 2. For d — "special" — 
2 cos j^, i a positive integer ( , )d is positive semidefinate. For other values 
o/cieR\0;(,)d has mixed signs. □ 

For d = 2 cos ^ define the radical Ra C TL^ hy < Rd, TL^ >d= 0. The 
radical Rd has first non-trivial intersection with the {£ + grade where 
it is 1— dimensional: Rd H TL^^^ = and Rd fl TL^^^+i = spanc(p^+i). The 
elements pi^i belonging to TL_^^_|_i (the £ + 1 grade of the generic algebra) 
arc called the Jones- Wenzl [W] projectors and a simple recursive formula for 
these is known. 

In this paper we will be particularly concerned with p^ and p4 (for d — y/2 
and ^-y^ respectively) which can be computed (from the formula on page 18 
[KL]). 
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, Pi generates a proper radical for d = \,— \; 

n 



u, u 



d U 



^ -1 \n ^ ( -1 \n 

genrerates aproprer radical for d = ±-v/2 , and d — 0; 



d'-2 



-d'+l /U 



d'-2d \n 



u ^ r^^u u^y 
n; d'-2d In Vi "^(Y n 



uu d KijJ , u u\ I 1 Vuy 
d^-sd'+i nn c/^-3rf^+2 \n n VTh/ d'-sd'+i 



A genrerates a proper radical for = ± 



'i+Vs' 



Figure 2.3 

It is known that when d is special, the ideal J(p^+i) generated by p^+i in 
TLd (the specialized TL algebra) is Rd- The notion of ideal closure in different 
algebraic contexts is essential to all that follows so we will be exphcit here; 
J{pi+i) is the smallest subset of TL^ containing p^+i so that if Ci, C2 G C;a,b 
belong to the subset and x,y E TL^ then: CiO + C26, ax, xa, a®y, and y ®a 
all belong to the subset. J(pf+i) is linear subspace of TL^ and is a two 
sided ideal under . and ®. So J is, by definition, closed under formal linear 
combination and all types of picture stacking: top, bottom, right, and left. 

For d special the algebra TL,^ contains many other ideals besides Rd (e.g. 
the ideal generated by diagrams with at least two "horizontal" arcs) but we 
find that when we move to the category, R^ becomes unique (sec Appendix). 
This motivates the definition of the Temperley-Lieb category TL^. 

The generic Temperley-Lieb category TL'^ is a strict monoidal tensor cat- 
egory over C(A) with objects A^o = {0, 1,2,.. .} thought of as that number 
of marked points in the interior of a horizontal interval. The indetermi- 
nate A determines d, above, by the formula d = —A^ — A~'^. The mor- 
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phisms Hom(m, n) is a C{A) vector space spanned by all pairing of the 
n + m points that can be realized by disjointly imbedded arcs in a rect- 
angle for which the m points are on the top and the n points on the bottom 
edge. The only difference from the algebra is that we do not demand that 
a nontrivial morphism have m = n. Again composition (") is vertical stack- 
ing and ® is horizontal stacking. The involution, the specialization of d 
and the notions of "ideal" are defined using exactly the same words as be- 
fore. Now the Markov trace < a,b >= tr(a6) becomes a Hermitian pairing 
Hom(m, n) x Hom(m, n) — > C. Theorem 2.1 continues to hold with TL'^ 
and TL^ replacing TL and TL^ respectively and for d special the radical Rd 
is still the ideal closure of pe+i- But in the categorical setting there is a new 
result conjecture by the author and proved by Goodman and Wenzl (see Thm 
3.3. in the appendix to this paper), which when combined with Theorem 2.1 
yields. 

Theorem 2.2. (Goodman, Wenzl) For d = special value = 2 cos TL^ 
has a unique non-zero, proper ideal = Rd = J{pe+i) and on the quotient 
TL'^ / Rd the pairing (, )d becomes positive definite, li d ^ special value but 
is of the form d = a + a, a a. root of unity, a 7^ ±1 or ± i, then TL^ has a 
unique non-zero proper ideal, the pairing (, )d descends to the quotient but 
has mixed sign. For other values of (i G M \ 0, TL^ has no non-zero proper 
ideal. □ 

We can continue to make the algebraic structure more flexible, more 
suited to both topology and physics, while retaining the key notion of "ideal" 
and uniqueness property set out in the proceeding theorem. One step in this 
direction is Jones theory of "planar algebras" [J2]. These are generalized 
categories with an operad structure replacing the notion of morphism. The 
TL-planar algebra, TL^ or TL^, if d is specialized, begins with a Hilbert 
space associated to an even number 2k of points marked on a circle: 
— span(imbeddable arc pairings in a disk D with 2k marked points on 
dD). To a disk with j— internal punctures D" and a relatively imbedded 1— 
manifold 7 C , where 7 has 2fcj endpoints on the A;*^^ interval boundary 

component and 2k endpoints on the outer boundary component, Jones asso- 

j 

dated (in an obvious way|^!) a homo morphism (^/i2fc, — ^ h2k- In the planar 

i=l 

algebra context the distinction between times (.) and tensor ((^) has been 
■^Let closed loops in D be assigned the multiplicative factor d 
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lost because there is no up, down, right, left. Instead we have "subpictures" 
of "pictures", i.e. restrictions of imbedded 1— manifold on a surface to a 
subsurface. 

Definition 2.3. A picture '-j in Y is an imbedded 1—suhmanifold (multi 
curve), proper ifdY A formal picture is a linear combination of pictures 
with identical boundary if 0. 

In Jones' theory there is no action by Dehn twist because surfaces are 
considered up to homeomorphism. 

We take a further step, and allow surfaces with genus > 0, here Dehn 
twist becomes crucially important. Consider an oriented compact surface Y, 
and the possible imbedded 1— manifolds ("multi-curves") 7 in Y. Picking a 
special value d for closed circles which bound a disk ("trivial circles") defines 
(i— isotopy. In section 1, we have defined ETL^^CK) to be the C— vector space 
of d— isotopy classes of closed null-bounding 1— manifolds modulo (i— isotopy, 
on a surface Y. 

Definition 2.4. Suppose a = T,aiji is a formal picture in a disk 6 C interior 
{Y) with fixed endpoints dji C dS. The ideal J (a) or (a) C ETLdiY) gener- 
ated by a are the d— isotopy of formal pictures of the form ax, x — TiXjXj, Xi 
a picture inY\5 with dxj — d^i, for all i and j, ax — EajXj(7iUXj). Dually, 

(a)* C ETL^{Y) are the functions annihilating (a). Concretely, y e (a)* iff 
y{ax) = for all x as above. The definition of ideal is the same in TLd{Y) 
and similar in the combinatorial settings: ETL^ {Y) and ETL'^ (Y) . 

One finds that the quotient ETL^/ J(p£+i) =: QE^ has (or better "recovers") 
the structure of a TQFT, or more precisely, a 2+1— dimensional unitary topo- 
logical modular functor (UTMF). For this, we must extend the definition of 

QEl to the case of a surface with labelled boundary {Y, t ). The essential 
feature is that QEl may be calculated by "gluing rules" applied to these 
smaller pieces. When we wish to emphasize the modular (UTMF) structure 
on QE^ we use the notation DEI = QE£ to recall the doubled 5*0(3) or 
"even" theory discussed in the introduction. Up the global |+) ^ |— ) invo- 
lution ~ on configurations, DE£ will be our model for the perturbed ground 
state space Ge,e, DEi = Gj^. 

A UTMF is a very natural way to model the topological properties of a 
two dimensional particle system without low lying modes in the bulk. Know- 
ing that the ground state has the structure of a particular modular functor 
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(DEi), tells us all the topological information about, excitation types, braid- 
ing rules (nonabelian Barry phase), 6^— symbol, 5"— matrix, and fusion rules. 
It is this structure that we have been seeking. 

The statement DE£ = QE£ is a purely topological one and it is possible 
to piece it together from the topological hterature using [BHMV] , [Prz] and 
[KL] . An exposition [FNWW] of the easiest modular functors is in progress 
and will explicate this isomorphism and contain a proof of theorem 2.5 below. 

But let us take a step back and explain this structure (UTMF) in a context 
where the gluing rules are obvious. Then we will summarize the axioms and 
finally explain the labels, pairing, and cutting/gluing operations in DEi in 
terms of functions on pictures. 

Let M(Y) be the vector space spanned 1— submanifolds (= pictures) 7 
of Y with no equivalence relation. Suppose Y is cut into two pieces by a 
circle a <Z Y , Y = YiUa Y2. The uncountable set X of all finite subsets of a 
will be the "labels" or "superselection sectors" of this theory. Neglecting the 
measure zero event that 7 and a are not transverse, we can formally write: 



where M{Yi,x) is the vector space of 1— manifold in 1^ meeting a = dYi, in 
the finite point set x. Equation (2.6) is the essential feature of a TMF as 
used by Wit ten [Wi] and formalized by Segal [S], Atiyah [A], Walker [W], and 
Turacv [T] . Many enormous "classical spaces" have this kind of formal struc- 
ture but it requires beautiful algebraic "accidents" to find finite dimensional 
"quantizations" of these. 

Bounding or "even" pictures span another (huge) vector space EM(Y). 
Let us set d = 2 cos and constrain the functions, in M*{Y) and EM*(Y), 
first by the d— isotopy relation and then annihilation by the ideal gener- 
ated by Jones- Wenzl relation pe+i. This yields the following quotients and 
inclusions in lines (2.2) and (2.3) 



M{Y) = ^M{Y,,x)(^M{Y2,x) 



(2.1) 



DK(A = ie'^^/2^+^) 
DK* {A = ie"^/2^+^) 



- TL,{Y) ^ M{Y), 
{p,+,Y^TLl{Y)^M*{Y) 



(2.2) 



DEI i — ETLdiY) i — EM{Y) 
DEt = {pi+iY ^ ETL*a{Y) ^ EM\Y) 



(2.3) 



25 



Theorem 2.5. The annihilating subspace (pe+i)* of M*(Y), we wrote it 
DK* , is in fact, the Drinfeld double [Dr] of the unitary topological modular 
functor (UTMF) derived from the Kauffman bracket at A = ie^'^l'^^'^'^ . This 
is true even at odd levels, I = odd, where the undoubled Kauffman bracket 
MF is flawed by having a singular S— matrix. DEi is a UTMF for i 2 
mod 4 and in these cases is a trivial double: V* ®V . 

Remark 2.6. For the corresponding even space, DEi the same MF arises 
for A, iA, —A, and —iA; A = the notation agrees with the intro- 

duction. 

Remark 2.7. The Kauffman bracket TMF (or TQFT), constructed in 
[BHMV], is not identical to the TMF derived from SU{2). In physics there is 
the loop group L{SU{2)) approach and in representation theory there is the 
quantum group {U sl2,q) approach and these lead to the same representation 
categories. Globalization of these representation categories (this view point 
is explained in [Ku]) yields the same MF. The pictures underlying the Kauff- 
man bracket are unoriented arcs. The Rumer- Teller- Weyl theorem shows 
these almost correspond to Rep q{SU{2)). However an important minus 
sign, the Frobenius-Schur indicator, corresponding to the quaternionic (not 
real) structure of the fundamental representation is missing. This minus sign 
propagates into the S matrix making the K and the SU (2) (TMFs) distinct. 
A different microscopic model which allowed arbitrary 1— manifolds (not just 
bounding 1— manifolds) could, depending on the local details, lead to VKi 
or DSU {2)£ so solid state physicists looking for anyons will need to be aware 
of this distinction in detail [FNWW]. The present models Ho/ and ad- 
dress only bounding 1— manifolds which correspond to (endomorphism of) 
the even symmetric powers of the fundamental representation which are all 
real. Thus Ki restricted to even labels EKi ^ S0{3)i, the SU{2) theory at 
even labels. The same holds, of course, for the doubles of these TMFs. 

Addendum 2.8. In [FLW2] it is shown that the braid representation of the 
"Fibanocci category"^ (F) is universal for quantum computation. DE3{A = 
e'fVio^ has 4 labels 0, 0; 0, 2; 2, 0; and 2, 2 and is isomorphic to F^F* implying 
that DE3 is also universal. 

Addendum 2.9. If Y has a fixed triangulation A then combinatorial ver- 
sions of the six vector spaces connected by maps (2.2) and (2.3) in Theorem 

^Greg Kupperberg's term for the even label sub-theory of {SU(2),3) also called the 
S'0(3)-theory at level 3. 
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2.5 are defined. Provided A is sufficiently fine, no information is lost; the 
left most combinatorial spaces are actually isomorphic to DKl and DEI re- 
spectively. The proof is the same as in the main topological theorem of [Fl]. 
Furthermore an estimate on the required fineness of A (it is linear in t) can 
be extracted from that proof. Also if Y has boundary and labels t (see the 
discussion of labels which follows immediately) are specified, then the left 
most combinatorial spaces are again defined and these map onto the TMFs 
with the given boundary labels t . 

To appreciate the last statement we set out Walker's axioms [Wa] for a 
UTMF. Fortunately these can be abbreviate due to two simplification: 1) 
the theories are unitary and 2) both are quantum doubles (i.e. the endo- 
morphisms, of another more primitive UTMF), so the central charge c = 
{c{V®V*) = c{V) + c{V*) = c{V)-c{V) = 0. Thus no "extended structures" 
or projective representations need be mentioned. For a concrete appreciation 
of these examples, see figures 2.4 and 2.5 where the particle types (= labels), 
fusion algebra, and braiding, and S— matrices in the cases DE3 are given. 

A labelled surface F is a compact oriented surface Y possibly with bound- 
ary, each boundary component has a base point marked and a label teC from 
of finite label set C with involution ^ containing a distinguished trivial ele- 
ment 0, fixed by^. For Kauffman SU{2), and S0{3) theories the labels are 
self dual a = a but we include the hats in the formulas anyway. A UTMF 
will be a functor V from the category of label surfaces, and isotopy classes of 
diffeomorphisms (preserving labels and base points) to the category of finite 
dimensional Hilbert spaces over C and unitary maps. 

Axiom 1 (disjoint union): V{YiUY2, tiU^a) = V{Yi, ti, ) <S)0^2: ^2), the equal- 
ity is compatible with the mapping class groupoids: 

V{hUh) = V{h)(^V{h). 

Axiom 2 (gluing): If Yg is obtained by gluing Y along dually labeled {x,x) 
boundary components then: 

V{Yg,t)^ V{Y,t,x,x) 

{x,x) € labels on the 
paired circles 

The identification is also compatible with mapping class groupoids - and is 
associative (independent of order of gluings). 
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Axiom 3 (duality): V{Y,t) = V{—Y,t)*, where — is orientation reverse on Y 
and ^ on labels, and * denotes the space of complex linear functionals. The 
Hermitian structures on V give vertical maps and the diagram below must 
commute: 



V{Y) i — >V{-Y)* 

I t 
V(Y)*< — >V{-Y) 

All these identifications are compatible with the mapping class groupoids. 
The Hilbert space pairings are compatible with maps: 

{x,y) = {V{f)x,V{f)y), where x,yeV{Y,t), 

and disjoint union: («i <^2, A /32) >= {<^i, Pi){o:2, ^2)- Writing 
a, j3eV{Yg, t) as a = and /? = according to axiom 2, then {a, 13) = 

X X 

^Sx{oix-i Px)- Where Sx — W Soxi- The symbols 5*00;; are the val- 

ues of a fixed function £ — > C \ {0} which is a part of the definition of V. 
Experts will recognize So^xi as the (0, i) - entry of the S— matrix of V: this 
is the map that describes exchange of meridian and longitude of a torus in 
the natural "label" bases. 
Axiom 4 (Empty surface): \/(0) = C 

'C, a = 



0, ay^O 




Axiom 5 (Disk): Let D he a disk, V{D, a) = 

Axiom 6 (Annulus): Let A denote an annulus. 
Then 

V{A,a,b) ^ 

To complete these axioms to a theory incorporating 3— manifolds Walker 
adds axioms 7 — 10. We will not need these here except to note that a 
3— manifold X determines a vector Z[X) belonging to V{dX). If X = F x /, 
V{X) = id e Hom(V(r X +1), V{Y x -1)) = V{Y) (g) V{Yy =: V {V{Y)). 

In the SU (2) theories it has been known since [Wi] that if X contains a 
labeled link "Wilson loop" (or suitable labeled "trivalent graph") then this 
pair also defines an element of T>V{Y). The simple idea is to regard the 
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1— manifold 7CF = FxOasa link labeled by "1", the 2— dimensional 
representation of s\{2,q) inside X = Y x [—1,+!]. This defines a map 
M{Y) — > Di{Y). If 7 is null bounding (in Z2— homology) on Y then there 
is a subsurface Yq C Y with d(Yo) = 7. Let G be a generic spine (trivalent 
graph) for Yo. Derived from Witten's theory and it bracket variation are 
combinatorial recoupling rules {6j symbols) which are exposited in detail by 
Kauffman and Lins in [KL]. We have adopted their notations (which caused 
us to rename Walker's trivial label "1" by "0") except in the choice of A, 

= q = e^'^*/^+^. To make d positive we choose A = ie'^*/^^"'"^, that is i 
times the primitive 4{i + root chosen in [KL]. For i = even our A is still 
a primitive 4{i + 2)*^^ root of unity, for i = odd, it is a primitive 2{i + 2)*^ 
root of unity but still defines a nonsingular TUMF on the even labels. 

Applying recoupling, 7 yields a formal labelling of G in which only even 
labels - odd dimensional representations — appear. This means that the set 
of possible morphism Z{X, G) is isomorphic to the endomorphism algebra of 
the sum of even labelled blocks. Restricted to even levels, the 4 choices for 
A differing by powers of i give the same Kauffman bracket UTMF and this 
agrees with the ^0(3) UTMF, which we call DEi. 

The recoupling relations on labelled trivalent graphs G G Y x 0, i.e. the 
6j symbols, are consequences of projector the relation pe+i applied to formal 
1— manifolds (and conversely p^+i follows from 6j). What is less direct [Prz] 
is that on a surface Y, pe+i alone generates the same relation as including 
y X C y X [— 1,+1] and then employing both p^+i and the Kauffman 
bracket relation X y4)(+A~^^. 

Abstractly we know the label sets for VKi and DEi, but we need to 
interpret these labels in VKi := TL^/ (pt+i) and DEi := ETL^/ {pi+i) resp. 
and in this context of recover the gluing formula. From a physical point of 
view it would be surprising if we could not localize because we expect the 
Hamiltonian H^ £ to define a stable topological phase for which the supers- 
election sectors of excitations define the label set. But such reasoning is in 
the end circular; it is better to have a mathematical proof that the candi- 
date ground state space G^/ has the structure of a UTMF and view this as 
evidence for or a "consistency check" on the physical stability of G^/- 

We now explain the "labels" for the theories DKi and DEi in terms of 
"pictures". A conceptual point is that the label has a kind of symplectic 
character: "half the label's information is a non negative integer < i which 
counts "essential" strands of 7 passing inward from a component C C dY. 
Think of this as "position" information. (Any "excess" strands correspond 
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to a descendent field or gapless boundary excitation.) The other half of the 
information ("momentum") is expressed as a symmetry condition on 7 in 
the bulk Y. The formal picture 7 must lie in the image of certain minimal 
idempotents — certain eigenspaces or projector images as constructed below. 

Abstractly, the label set C for DEi = QEi := ETLd/{pe+i) may be 
written as: 



£ = {(0,0); (0,2); (2,0); 



+ 1 



+ 1 



} 



The "position" part of the doubled label t = {a, b) for 7 G QE£ on a bound- 
ary components C C dY is |a — 6|. This quantity is the smallest number ^ of 
domain wall (7) intersections with C , ^ — \a — b\, as C varies over all imbed- 
ded loops parallel to C (i.e. cobounding an annulus with C) and the domain 
wall 7 also varies over all (p^+i)— equivalent pictures. The "momentum" 
part of the label is an eigenvalue. 

Let us do this more carefully. We follow [BHMV] to define what Walker 
calls an "annulus category" Af. A — S^xl, obj(A^) = {set of even number 
of points on S^} then an clement of morph (A^) are all formal combinations 
of pictures in A, which beginning on the object in S*^ x — 1 and end on the 
object in 5*^ x +1, and which obey the relations: rf— isotopy and pe+i. (Recall 
(pi+i) = negligible morphisms of TL^. Also see appendix.) Suppose that Y 
is a surface with connected boundary dY — C, then there is a gluing action 
of Af on DEi{Y) : / e morph (a^), g G DEi{Y), f og e DEi(Y); 
f o g[xi (E) Zj) = Xi^j) where the coefficient of / on the picture Xj is Xii 
fi^i) = Xi ^iid g{zj) = ujj. For this action to be defined we must pick 
an identification Y = Y U A. Also C has a fixed parameterization and 

an orientation; these tell us which end of A, x —1 or x to glue 
to C. Technically, this means one of Af "''^ or Af is acting according to 
orientation. Since we will not make calculations, we will not be careful in 
choosing orientations and in distinguishing categories and their opposites. If 

Y has k boundary components the A;— fold product acts on DEi{Y). 

k 

The reader should note that in the cases where there is a miss match between 
boundary conditions on dY and there is by definition no action defined. 

k 

This would not make sense if we were dealing with algebras and is precise by 
the extra fiexibility that make linear categories, sometimes called algebroids, 
a useful generalization of algebras. 
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So Aa/ =: C has an action or — to use the usual terminology when 

actions are linear — a representation on DE£{Y) —: V. 

To clarify, for each 9— condition on ]JA there are two corresponding finite 

k 

dimensional vector spaces DE£{Yi^to a a-condition) and 
DEiiXont A 9-condition) =: Kut- A morphism in7out e morph(C) induces (by 
gluing) a linear map : \4n — > Vont- The construction is so natural that all 
the required diagrams commute, and gluing, indeed, defines a representation. 
(There are some technical points but these are well understood and will cause 
us no trouble. As annular collars are added to Y various "association" must 
be chosen so the action is "weak" not "strict" . Also it is sometimes conve- 
nient to forget the parameterization of C and only remember the base point 
* = 1 e S*^, and further to replace the uncountable object set — finite sub- 
sets of S^— with the countable object set consisting of one exemplar object 

for each finite cardinal 0, 1, 2, This processes is called skeletonizing the 

category.). 

The positivity properties of the pairing ( , ), figure 2.2 and Thm 2.1, and 
its extension implies that its finite dimensional representations decompose 
uniquely into a direct sum of irreducibles. The arguments for this are nearly 
word for word what is said to prove that a finite dimensional C*— algebras is 
isomorphic to a direct sum of matrix algebras. The algebroid context changes 
httle. 

So let us decompose V as a representation of C into a direct sum of its 
irreducibles. We record multiplicities by tensor product with a vector space 
Wi on which no action of C exist: 

irreps. C 

The index i is a multi index i — {ii, . . . ,ik) and counts the "admissible la- 
beUings" of dY — C-JJ . . .UCk. In the theory DEI, the possible components 
of i are the 

{irreps of } =: £, (2.5) 

the label set of DEI. The involution ^ is induced by orientation reversal on 
A and conjugates representations: it happens to be trivial in these theories. 

A final categorical comment: the form of r.h.s. 2.4 suggest it represents 
a 2— vector: a linear combination (in this case of irreps.) with vector space, 
rather than scalar coefficients. This is, in fact, the correct categorical setting 
for MF{Y) when Y has boundary. 
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So far this discussion of the label set has been rather abstract but it is 
possible to make explicit calculations by considering the annular categories 
as operators on TL^, the Temperley-Lieb, or "rectangle" categories. A label 
a e irrep (a^) is generated some idempotent a = k^k = morph (k, k) C Af, 
which is a linear combination of annular pictures. Quoting a result which will 
appear in [FNWW] with full details, we describe (Figure 2.4) the idempotent 
a for the four labels of DE3, a = (0, 0), (0, 2), (2, 0) and (2, 2). In the language 
of rational conformal field theory these labels are the 4 primary fields . They 
are given below as formal pictures in annuli. 

Previously, we only considered ideals to be generated by formal pictures in 
a disk. But now we can let a, by stacking formal pictures in annuli, generate 
an "annular ideal", J. Elements b E J may have more than k boundary 
points on 5*^ x ±1; such b are the descendent fields. 
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0,0 



2,2 



where D - 



m = - 




in the quantum dimension of DEli and 



'd 



U 

n' 



d =-A''A-\A = ie'""\ 



, 1 2 . 37r 




where uj^ = [(), j , = (^2, 2j , and the crossing symbol on 
strands labelled ^ is expanding the Kauffman bracket; 

^Yd n^"/^^+/^'nn 



2, Oj is the complex conjugate of (^0. 2 j . 

Figure 2.4 

The idempotent (0, 2) has among its many terms five principal terms, 
pictures containing only arcs going between inner and outer boundaries of A 
and no arcs which are boundary parallel. The other terms enforce orthogo- 
nality of (0,2), to the descendents of (0,0) and (2,2). The principal terms 
written out below. 



(0, 2) :: / + e^^^'^F + e^^'^^F^ + e^'^/^F^ + e^^^'^^F^ + lower terms F (2.6) 
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fractional Dehn Twist: F =©, / = 0. (0, 2) is a e^'^*/^ eigenvector of F. 
The powers of F are obtained by radial stacking of annuli. 

In the case [a, — 6j| = 2fcj, /cj 7^ 0, consider the commuting actions of 
F = e'"*/*^' twists of Ci on QE^. The resulting eigenvalues turn out to be 
distinct within the "position" = \a — b\ summand ofV(Y). These eigenvalues 
add the "momentum" information which determines the labels: the minimal 
projections to eigenspaces. 

For ki = 0, Dehn twist acts by the identity, so here prescription must 
be different. Suppose s is a configuration on Y which has constant spin 
"monochromatic" either |+) or |— ), near Ci and let [s] be its image in QEi. 
Define an action on [s] by adding an annular ring of the opposite spin in 

interior (Y) immediately parallel to Ci and define /5„ = E 2xR^ where 

x=0 ' 

consists of x parallel annular rings of the opposite spin stacked up paral- 
lel to Ci but in interior (Y), and where 5'^,^; is the S'-matrix = -^=sin|^ 

of the undoubled theory. For n = 0, . . . , [^J, the maps Pn are commut- 
ing projectors (up to a scalar), which also commute with twists around other 
boundary components. For fcj = "position" is refined to a label by applying 
the idempotent Pp. The trivial label is Pq. The image of Pp turns out to be the 
_j^^2p+2_|_^-2p-2-j eigenspace of the actions of R on the \a — b\ = "position" 

summand of V{Y). Thus when dY 7^ 0, DEi{Y, t ) is an orthogonal sum- 
mand of QEi{Y) determined by specifying a minimal even number rii, of arcs 
reaching each boundary component Ci (but not parallel to it), nje{0, ...,£}, 
and an eigenvalue Aj at Ci C dY. Note that gluing different eigenspaces im- 
ages automatically implies a trivial result as required by the gluing axioms. 
This is immediate from the commutivity of the minimal idempotents. For 

surfaces with boundary, we may write DEi{Y) := 0^yt)y 

admissible labelings t 

then DEi{Y) = QEi{Y) in this case as well. 

For the first computationally universal case, level i = 3, the S'-matrix, 
F— matrix (= 6^— symbol), the action of Dehn twist, and Verlinda (fusion) 
relations are listed below. The only interesting F— matrix in the undoubled 
theory occurs when all four external label = 2, so F reduces to a 2— tensor. 
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1= 3 

S - matrix = tensor square of: 










2 





V r sin'% 




2 




-Xfs ^^"/^ 


F- 


matrix = tensor square of : 






2 







2 


-d- d-y 
d-y d- 


^ i+Vs 

, d = 

2 



In the undoubled thoery the nonzero fusion coneiEcients are: 

^000 ' ^222 ' ^022 ' ^202 ' ^220 ' ^ ^he doubled theory, 25 
nonzero A' 's occur: N......i^i^. obtained by inter leaving the sub- 
scripts of any two A'^^.j and A';, J, above. The action of Dehn 

twist: , .,- 

Dehn (0,0) = (0,0), Dehn (0,2) = e'"" (0,2), 

Dehn (2,0) = e"" (2,0), Dehn (2,2) =(2,2) 

Figure 2.5 

We have come to a point where we can study the difference between span 
{2 — colorings (|+), |— ))} = Goi and the ~ — invariant combinations G^^. 
We may begin with the enhancement of ETL^iY) to ETL[{Y), manifold 
2— coloring modulo c?— isotopy (for d = 2cos7r/£ + 2). The enhancement 
leads the "color reversal particles" of figure 0.2 which do not fit exactly 
into the UTMF— TQFT formalism, (but perhaps a Z2— graded version?) as 
they do not raise the ground state degeneracy on the torus. They should, 
however, arise physically and contribute to specific heat. We will return to 
these shortly. 

First, we show that there is only one lifting to the enhancement of the 
projector relation pe+i : for i odd let p^]^f^ and pj^f^ denote the extension to 
2— colorings of pe+i by the relation applied to a 2— coloring in a neighborhood 
of a transverse arc which crosses i + 1 strands of 7 from black to black and 
white to white respectively. The reader may take black = |+) and white 
= |— ). If £ is even, noting that all the projectors Pi have left-right symmetry 
there is only one way to lift pi^i to ETL[. 
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Proposition 2.10. J {pf^^^) = J {pj+i^). 

Proof: We may use under crossings to indicate formal combinations of 
TL— diagrams which are consistent with the Kauffman relation: 

X-)(--n 

Figure 2.7 

Now consider (for £ = 3) the following sequence. 




Figure 2.8 



These 6 steps effect across the arc a with an apphcation of p^f^. 

□ 

For a closed surface Y let QEl(Y) be the enhancement ETL[(Y)/ {pY+i^) 
= ETL[{Y)/ (pj^f^) and QEl{Y) — > QEe{Y) the forgetful map. We do not 
know if this map is always an isomorphism. However for the case of most 
interest £ = 3, Proposition 2.11, below shows that QE\{T'^) ^ QE^iT^), 
the 2— torus. Since dim(K(T^)) = |£| the cardinality of the label set, this 
implies, that largest quotient of QE^ having the structure of a UTMF is 
isomorphic to QE^. 

Proposition 2.11. QE'^iT^) ^ Q3{T^). 

Proof: Let B be the black and W the white coloring of while M and 
are the one and two meridional ring colorings respectively: 
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B W M 




L D 

Figure 2.9 

Applying p4 across the arc 7 (and a short calculation using Figure 2.3) 
yields: = 3M — W . Since M and are black-white symmetric the 
third term must be symmetric as well, hence W = B. 

It follows quickly that QEl{T^) = C — span [W, M, L, D) where L is a 
longitudinal and D a diagonal ring. Forgetting the 2— coloring and retaining 
only the domain wall we get a basis for QE^{T'^) = DE3{T'^). □ 

It is possible to brake the color symmetry ~ by adjusting the Hamiltonian 
to fix the color = |— ) at some plaquet on each component of Y. This ad- 
justment creates a new ground state G canonically isomorphic to the former 
G~^, so we drop the from the notation. However this does not obviate the 
need to study the enhancement. The point is that localized color- reversing 
excitations remain and are expected physically. These, when realized on an 
annulus algebra, have opposite coloring on x —1 and x +1, and so 
cannot be glued into a ground state on T^. 

Let us see how this works in the simplest example, the level = 1 theory 
i = 1, d = 1, A = ie^'^''-!^'^ . When we make no "evenness" restriction this 
theory, Dl, is also called Z2— gauge theory [SF] and [Kl]. It has four labels: 
= 1(0 - i?), m = i(0 - R), e = i(J + T), and em = \{1 - T) where the 
pictures in these combinations are: 




Notice that the labels are orthogonal under stacking annuli. One may 
check that braiding m around e or em introduces a phase factor = —1, as does 
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braiding e around rn or ern. The even theory DE1{A = ie^'^'^^^^), has only one 
particle (0, 0)— which is the trivial particle, and has dimension = 1 on T^, and 
so is quite trivial. In quantum systems with other microscopies (e.g. [Kl]) 
can easily realize D{1) but in our set up the pictures do not arise directly but 
indirectly as a domain wall so / and T make no sense. However R does make 
sense as a domain wall between |+) and |— ) boundary conditions at opposite 
ends of A. In fact, we may define an elementary excitation of DE{1) at a 
plaquet c by using the local projector at (|+)c + \~)c) ® ((+|c + (— ||c), 
instead of the ground state projector (|+)c— |— )c) ® ((+|c — (~|c)- Thus 
the "m" particle can arise as an excitation of DEI even though it does 
not contribute to ground state degeneracy. This is the prototypical color- 
reversing particle. 

Regarding the other information in figure 0.2, the coloring preserving or 
"label" excitations (the irreps. of A/) are counted by the (even, even) lattice 
points in [0,£ +1] x [0,£ + 1]. The color reversing labels are (odd, odd) 
C [0,^ + 1][0,^+ 1]. The -S-matrix of the (undoubted) SU{2) or Kauffman 
theory when restricted to even labels is singular (precisely) for £ = 2 mod 4, 
for example at £ = 2 





1 V2 


1 
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V2 
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1 
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1 1 
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2 2 

1 1 

2 2 





When S is nonsingular, £ 1, 2, 4 and the number of braid stands > 5 it is 
known [FLW2] that the braid representations are dense in the corresponding 
special unitary groups. 

3 Perturbation and deformation of Ho/ 

As remarked near the end of Section 0, excited states, i.e. anyons, will be 
studied as ground states on a punctured surface with labelled boundary. In 
the large separation limit, the braiding of anyons can be formulated as an 
adiabatic evolution of the ground state space on a labelled surface Y with 
boundary. So in the present section we confine the discussion to ground 
states. Although boundary is assumed to be present and labelled we will 
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nevertheless consider only perturbations acting in the bulk so the role of the 
boundary is peripheral in this section. 

The passage from Go/, the ground state space of Ho/, to the deformed 
ground state space G^/ of H^^ does not result from the breaking of a sym- 
metry, in fact Go/ has no obvious symmetry. Rather it is the creation of 
new "symmetry" : topological order. If a perturbation V is breaking an ex- 
isting symmetry then only the original ground state and the effective action 
of V at the lowest nontrivial order is relevant. But in the present case, to 
understand the effect of a perturbation eV , one should first describe all low 
lying (gapless) excitations above Go/ and then see how V can act effectively 
on Go/ through virtual excitations. For example in the toric codes [Kl] the 
ground state space may be rotated in an interesting way if a virtual pair of e 
(electric) particles appear, tunnel around an essential loop (of combinatorial 
length = L), and then annihilate. In the case of toric codes there is an energy 
gap to creation of (e, e) pairs so the above process has exponentially small 
amplitude in the refinement scale L ~ e~^/^° . In contrast, for level £ < 2 
we expect the ground state space Go/ to be gaplessQ(in the thermodynamic 
limit) and processes which act through virtual excitations will be important 
in the perturbation theory because excitations are cheap. However it seems 
hopeless to analytically describe these gapless excitations so we skip this step 
and resort to an ansatz (3.4) stated below. It asserts that G^/ is modelled 
as the common null space of local projectors acting on Go/- We argue for 
this via an analogy to FQHE, uniqueness considerations and "consistency 
checks" . 

From section 2, the reader knows that we wish to identify G^/ with 
Go// (pi+i) (for suitable values of e), and this is what the ansatz implies. 
The fact that Go//{pe+i) = DE^ (see § 2 and [FNWW]) has the structure of 
an anyonic system (mathematically a UTMF) is the first consistency check. 
There will be one more presented in section 4. Let us prepare to state ansatz 
3.4 carefully. 

Definition 3.1. An operator (9 on a tensor space Ti = (^C^ is k— local 

if it is a sum of operators Oi each acting on a bounded (< k) number of 
tensor factors and id on remaining factors. We say O is strongly local if the 
index set {V^} are vertices of a triangulation A and O = T,Oi, where each Ot 

i 

is fc— local with the k active vertices spanning a connected subgraph of 
^For £ — 3 the gap may be extremely small as explained latter in this section. 
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A. All Gi are assumed isomorphic and with fixed isomorphisms Gj — > Gj 
inducing isomorphisms Oi = Oj. In the latter case, we call {A,0) a quantum 
medium . 

Note 3.2. In the special case that a family of strongly local operators {Oi} 
are projection onto 1— dimensional subspaces, the system {Oi} is equivalent 
to what topologist call a combinatorial skein relation ([L] [KL]), though in 
the topological context equivalence of isotopic pictures is implicitly assumed. 
An example of a (14 term) skein relation is p^+i = (see fig. 2.3), apphed 
to 7, the dual— 1— cell domain wall between |+) and |— ). A skein relation is 
a local linear relation between degrees of freedom, domain walls in our case. 
The intersection of all the null spaces f]null((9j) = null((9) is the subspace 

i 

perpendicular to the equivalence classes in 7i defined by the combinatorial 
skein relation. 

Definition 3.3. The joint ground space (jgs) of {Oi} is H-^o.i where E^^i is 

i 

the eigenspace corresponding to the lowest eigenvalue Ao of Oi. 

The jgs {Oi} is not necessarily the lowest eigenspace of (9 = SOj because 

i 

jgs {Oi} can easily be {0}. In this case the Hamiltonian is "frustrated". It 
may happen that O has long wave length excitations at the bottom of its 
spectrum which do not show up in the spectrum of Oi. However if O defines 
a stable physical phase it is an optimistic but not unrealistic assumption 
that jgs {Oi} = ground state space (O). For example, this occurs in the "ice 
model" or "perfect matching problem" on the honeycomb lattice [CCK] and 
in the fractional quantum Hall effect (FQHE). 

The FQHE begins with a "raw" state space H, the lowest eigenspace for 
an individual electron confined to a 2— dimensional disk D and subjected to 
a transverse magnetic field B. This Ti. is called the lowest Landau level . Each 
level can hold a number of spin + and — electrons ~ area D/ (magnetic length)^ 
and the fraction of that number actually residing at the level is called the 
filling fraction u. In a spherical model, the Coulomb interaction H between 
pairs of electrons (cj, ej) can be written [RR] as a sum of projectors onto var- 
ious "joint angular momentum subspaces" pij{{2k + 1)N^). The null space, 
null(ifq), for 

g-3 
2 

H, = J2J2P^Ai'^k+l)N^) (3.1) 

fe=0 i<j 
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is nontrivial and, of course, is the joint ground space jgs of the individual 
projectors in the sum. In fact, nuU {Hg) is Laughhn's "odd denominator" 
state space at u = K 

Ansatz 3.4: For weU chosen e, the perturbed ground state space wiU 
be stable and can be written as = jgs {\si) {si\)nGo for some strongly local 
family of projectors acting on H. Equivalently if s° is the orthog- 

onal projection of Sj into Go and : Go — > Go is the corresponding 

projector then G^ = jgs (|s°)(s°|). 

In topological terms the ansatz asserts that the reduction 
Go — Ge occurs by imposing a skein relation. The ansatz is essentially 
a strong locality assumption. 

As discussed above, the Laughlin ^ = ^ states, q odd, follow this pat- 
tern with the Coulomb interaction between electrons playing the role of the 
perturbation on the disjoint union of single electron systems. Since the Lan- 
dau level has no low lying excitation the analogy is closest with Go/, £ > 3. 
Theorem 2.5 gives us the following: 

Implication of Ansatz 3.4. Suppose is subjected to a sufficiently small 
perturbation £ < 2 or an appropriate deformation, i >3, which partially lifts 
the log extensive degeneracy of the ground state Go,^ to yield a strictly less 
degenerate ground state G^/- If we assume the stability of G^/ we expect 
Ge/ to be modelled as Ge/ = Go//{pe+i) = the modular functor DEi. 

For the projector p^^i to arise as an effective action of eV on Go/ = 
C < d— isotopy classes of domain walls >, for e > 0, various sets of £ + 1 
walls must have a polynomially large probability of simultaneously visiting 
the support Ui of some local combinatorial Oi — p\^i enforcing pe+i. The 
walls must visit a site Ui or pl^^ cannot enforce orthogonality to the relation 
vector pi^i (as depicted in figure 2.3 for £ < 3.) 

The notion of a combinatorial instance of p^+i was developed in [F2] . It 
amounts to a discretization of the smooth domain wall diagrams (Figure 2.3) 
by choosing specific superpositions of local plaquet spin configurations (with 
the spin state of the plaques at the boundary of the configuration constant) 
to represent the smooth relation p^+i. Evidently there are many distinct 
combinatorial patterns which are instances of a fixed pe+i. The simplest of 
these amount to geometric rules for simphfying £+1 domain wall when these 
run parallel for (roughly) £ + 1 plaques. As discussed in [F2], imposition 
of such a combinatorial relation in the presence of mild assumptions on the 
triangulation A, produces a result isomorphic to the smooth quotient. It is 
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sufficient that the triangulation must have injectivity radius » £ + 1 and 
bounded valence. So A should subdivided (to approach a thermodynamic 
limit) as shown in Figure 3.1. 

Heuristic: There is a curious pattern observed in Figure 2.3 (and further 
computer calculation of Walker (private communication) , for ^ + 1 = even 
and d — 2 cos j^, the sum of the coefficients in Pi+i, in the geometric basis 
{g} is zero, and for i+ 1 = odd, the sum is small. The geometric pictures {g} 
may be filtered by an integer weight n which counts the fewest sign changes 
on plaquets — in topological terms, the fewest domain wall reconnections on 
"surgeries" — required to transform the straight (identity) picture to g. So, 
referring to Figure 2.3, the first term for has weight = 0, the second term 
has weight = 1. The terms for have weights 0, 2, 2, 1, 1 respectively. Notice 
that sign (coefficient (g)) = (-l)™ght (f). 

This suggests that V = E o"^ is a reasonable choice for V to obtain 

c, plaquct 

figure 0.1. The Pauli matrix has +1— eigenvector |+) + |— ) and —1— eigen- 
vector I +) — I — ) , and thus assigns a lower energy to combinations of geometric 
pictures g which have a (—1) phase shift associated to domain wall surgery. 
The perturbation 

V — S(|-|-) + |— ))((+! + (— I) = ^T,{id + a^) contains terms which annihilate 

c c 

antisymmetric combinations of approaching domain walls of 7: (|)() — 

Nonzero entries coupling all the terms of pi^i occur first at order i, i.e. 
in V^, so one may expect that this is the order at which an effective action 
arises on Go/. 

It is now time to treat the statistical physics of a general ground state 
vector ^ G Go/. A perturbed Hamiltonian H^/ will not have a ground 
state modelled by Go// (pe+i) if the domain walls of ^ have an effective 
tension. The simplest place to see this is on a closed surface Y with the 
triangulation A determining a metric. When the domain walls 7 C y are 
pulled tight under tension they will stand a bounded distance apart and have 
exponentially small amplitude for simultaneously entering a small locality Ui 
so p^^^ will be unable to act. 

Measuring any ^' e Go/ via the community family {cr^} projects ^ 
into the geometric basis. This results in a probabilistic spin configuration 
which is Gibbs with probabihties proportional to n# where n — d'^ — 
(2cos7r/£ + 2)2. 

Let us write ^ G 7i, |^| = 1, in the classical basis of spin configurations, 
^' = Safcl^'fc). We say, consistent with measurement of any observable which 
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is diagonal in the j^^) basis, that has probabihty la^p. Thus becomes 
a random classical component of the random configuration, meas.(\E'). Just 
as one asks about the typical Brownian path, we ask what a typical 
looks like. There will be a competition between energy and entropy. Since 
d > 1, the Hamiltonian Ho,e "likes" trivial circles and will place a high weight 
on configurations with most of the surface area of Y devoted to a "foam" 
of small circles. However entropy favors configurations with longer, fractal 
loops which exhibit more variations. From the critical behavior of loop gases 
we know that for d < \/2 entropy dominates and ^ G Go/ is a critical 
Gibbs state with typical loops fractal. For d > y/2 energy dominates and 
the Gibbs state is stable: To free up dual lattice bonds to build this foam 
the topologically essential part 7"*' of 7 will be pulled tight by an effective 
"string tension" . 

Recall from the introduction that the Gibbs weight on a loop gas state 
7 is proportional ^(7) = e"'^^*"*^^ '^"^th 7)^* components 7^ our is the self dual, 
k — case. 

Let us be explicit. In any ground state vector ^ e Go/, ^ = Sai\&i, we 
have seen that the coefficients ai and 02 of d isotopic configurations and 
^2 satisfy 01/02 = c?*Y^*^ where #i is the number of trivial domain wall 
components, " trivial loops", of ^j. The Pauli matrices cr^ = | q -1 1 apphed 
at vertex = v form a commuting family of observables so we may "observe" in 
the geometric basis of classical spin configurations to obtain meas.(^) 
and we see that the ratio of probabilities of observing ^1 verses ^2 is: 

2/ 

Thus observing ^ in the classical basis yields a Gibbs states 
meas.(^) :: e'/^^^"^^) for E{^i) = -#i and P = 21ogrf. 

Such probabilistic states are called "loop gases" and have been extensively 
studied [Ni], e.g. in the context of the 0(n)— model. It is beheved that, in 
the self dual case k — 0, d < V2, £ < 2, there is no string tension and that 
the Gibbs state is critical: typical loops are 1 /polynomial in size and cor- 
relations decay polynomially. Furthermore the familiar "space = imaginary 
time ansatz" (see lines 3.7-13) suggest that this regime, i < 2, should have 
Go/ gapless. For i > 2, the loop gas is beyond the critical range. For these 
values correlations of the loop gas decay exponentially and it is believed the 
loops are in "bubble phase" where any long loop forced by topology would 
be pulled tight by an effective string tension. The corresponding Ho/ i > 2 
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should be gapped, above its (polylog) extensively degenerate ground state 
space Go/ (compare with line 3.8). It is in this case, specifically for £ = 3, 
that we still may hope the ansatz describes G^/ for some family of deforma- 
tions H^/, ei > e > eo > as suggested by the phase diagram. Figure (0.1). 
This would imply G^/ = Go// < pe+i >= Go//Ri = DEg. 

The gap above Go, 3 and therefore eo might be quite small. A loop gas 
with k = (defined on a mid- lattice see Ch. 12 [B]) is closely related to 
the FK representation of the self-dual Potts model at g = = d'^. For 
i = 3, d = and q = 5.6. Although the self-dual lattice Potts model 

in 2— dimensions have second order transitions (critical) for g < 4 and first 
order transitions (finite correlation length) for g > 4, exact calculations show 
that for q = 5.6 the correlation length ( though finite is several hundred 
lattice spacing [BJ]0. 

Our "alternative" Hamiltonian H'^ ^ resulted from an effort to sharpen the 
relation between meas.(^') and the (FK) Potts model. I would like to thank 
Oded Schramm for helpful conversations on this relation. Recall that H'^ ^ is 
a Hamiltonian on Hilbert space of spin = 1/2 particles on the bonds of a sur- 
face triangulation, cellulation, or lattice. We will work locally and so ignore 
contributions of the global Euler characteristic x(^) to the formulas below. 
Also, we write " =" to mean that the equation holds up to a fixed extensive 
constant, like the total number of bonds. Recall from the introduction that 
we consider the union of |+) bonds and disjoint from this the union of |— ) — 
dual bonds. Our "loop gas" is on the mid lattice separating the from 
the |— )— clusters. Let E be the number of edges and E* the number 
of |— )— edges, C(C*)) the number of clusters (dual clusters) with the con- 
vention that isolated vertices (dual vertices) count as clusters (dual clusters). 
Let L be the number of loops in the loop gas. The Potts model parameter 
q (number of colors), it turn out, should be set as g = d'' = (2 cosvr/f + 2)^. 
Finally, < p < 1, denotes a probability. 

We have two basic equations. Every loop (in the plane) is the outermost 
boundary of either a cluster or dual cluster so: 

L = C^C* (3.3) 

Also there is an Euler relation: 

C*" = "C + E (3.4) 

thank Steve Kivelson for pointing out the existence and relevance of these calcula- 
tions. 
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Now wc can rc-cxpress the loop gas Gibbs weight uj in terms of clusters in 
the FK Potts model: 



a;(spin conf.) " = "(d^)^ = {d'f^'' = {d''f{dTl'' " = " 



Because of the ~ - symmetry between |+) and |— ) we expect p to be the 
self dual point for this value of q, and this we check this below: 

using (3.4) we get p'^^{l — p)"^'^ = q^, so 

P = (3.6) 

we have proved. 

Theorem 3.5. Observing with the family {a^ for all bonds b} maps any 
ground state vector ^ e G'^^ into a Gibbs state meas.(^) of the self dual 
Potts model for q — {2 cos j^Y- □ 

This justifies thinking of H'^ ^ as a "quantum Potts" model and contem- 
plating a diagram of relations: 



conf. field theory, 
5^2, level £ 



correlation functions 



infra-red limit 



TT 



{SU, e)UTMF 



double 



DEi 



zero modes 



Potts q = (2 COS ) 



Potts q = {2 cos ) 



+ 2^ 



V 4 observe , 

< d 



quantum Potts 
(2 COS- 



Figure 3.1 

Let us return to the heuristic connecting the spectral gap above Go,e and 
the statistical properties of a ground state vector ^ G Goi- Let A and B 
be two strongly local operators on a quantum medium {A,H). For example 
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measuring A might be cr^ at plaquet c. Assume, first, that there is a unique up 
to a phase ground state vector Q and a spectral gap = 5 above = 0. 

Then if we evolve in imaginary time: 

"imaginary time correlation" = {Qe^" Ae~^" BQ) = {QAe'^^BQ). (3.7) 

Writing {BQ) = (QSfi)Q + *B and (QA]* = + note that non-Q 

summands and decay under imaginary time evolution at a rate > e~^^. 
Thus applying e~*^ to either the bra or the ket in (7) gives: 

"imaginary time correlation" — > {{flAfl)fl{flBfl)fl) = (flAfl) (flBQ) 

(3.8) 

where the convergence ( — >) is exponential. Denote the spatial translations 
of A (by Af), the analogous operator to A acting near a site at distance i 
from the support the original of A. The simplest expectation is that spatial 
correlations should behave in the same way as imaginary time correlations, 

"space correlation" = {QAeBn) "^P°'^"'^*'^^y , (3.9) 

This statement is precise in a Lorentz invariant context but is expected also 
to hold in greater generality provided that there is some linkage between 
temporal and spatial scales. 

A code space G C 7i is an important generalization of a 1— dimensional 
subspacc (sec [G] for examples and discussion in other notations). 

Definition 3.6. We say G dTi is k—code , with k measuring the strength of 
the encryption, if for any strongly k— local operator Ok the composition: 

G J^n^n^G (3.10) 

must be multiplication by some scalar (perhaps zero). Uq is orthogonal pro- 
jection onto G. 

The importance of a code space is that it resists local perturbations. In 
[G] Gottesman states, in different language his, Thm 3: 

Theorem 3.7. Suppose £ C is a subspace of a Hilbert space and S 
is a linear space, called "errors" , of the operators HOM {7i, H) so that the 
composition below, for any E & £ is always multiplication by a scalar c{E): 

L TC ^ H — ^ C Then C constitutes a "code space" protected 

from errors in £. This means that there is a physical operator (composition 
of measurements and unitary transformations) which corrects errors coming 
from E & S. (H^ is orthogonal projection onto C) □ 
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We have already encountered a code space. The space DE£ = Go,i/ {pi+i) — C 
Ti. has the code property (and similarly for G'^ ^/ (pe+i)). From the theo- 
rem and the disk axiom (section 2) we see immediately that the dual space 
Go,i n {pe+i}* C 7-^ is a code space for operators (errors) supported in any 
fixed disk D gY. This is true long before the refinement limit, we need only 
a modest level of refinement before the combinatorial quotient exactly as- 
sumes the structures of a unitary topological modular functor (UTMF) [F2] . 
Then axiom 4 (section 2) says that V{D, 0) = C and V{D, a) = whenever 
a label a 7^ 0. Thus any operator supported on D must act as a scalar. 

Observe that even when a ground state space G is degenerate, if it is 
nevertheless a code space and also has a spectral gap of 5 > between it and 
the first excited state, the argument for the decay of spatial correlations 

{QoAeBrio) (3.11) 

remains valid for any flo^G: 

\n,B) = {noBn,)n, + oQi + 0Q2 + • • • + *s ± g, (3.12) 



{QoA\* = {noAno)n^ + or^i + 01^2 + • ■ ■ + ^a, ^ g, (3.13) 

{Qo, fii, 1^2 • • •} an orthonormal basis for G. Applying e~*^ to say ket (3.12) 
and then pairing with (3.13) observe that (floAfli) =0,i > 0, so: 

"imaginary time correlation" ^^P°"^"^'^"y ^ {QoAQo){'^oBQo) (3-14) 

Thus the usual heuristic: gap < > finite correlation length, gapless < > poly- 
nomial decay, is not less valid for code spaces than simple, non-degenerate 
ground states. Hence, the expectation is that code spaces G are protected 
by a spectral gap and meas.(^), ^ E G has finite (or even zero) correlation 
length. 

Curiously, it is the square of the Beraha numbers, n'^ = d"^ = q = 

(2cos^j , i > 1, which enter as the weight of a circle in the loop gas 
Gibbs state. This means that for £ > 3 these systems are outside the crit- 
ical range in the thermodynamic limit. But recall [BJ] that, for ^ = 3 the 
resulting theoretical stability is extremely weak. 

We would like to propose the possibility that the transition Uc from critical 
to bubble phase on a surface Y might be sensitive to roughening, i.e. an 
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increase of Hausdorff dimension. Roughening appears to increase the entropy 
of long domain walls giving them more dimensions to fluctuate in, so one 
might expect the energy/entropy balance point to increase to ric > 2. There 
are two arguments for Uc = 2. One is an explicit study of the spectral gap 
of the corresponding transfer matrix in the 6— vertex model. This needs a 
geometric product structure and so will not apply to a typical rough surface. 
The second (ch 12 [B]) is topological and seemingly does apply. It used the 
Eluer relation to create a translation between the Potts model, a loop model, 
and an ice-type model. Formally, when the Potts parameter q crosses 4, 
9 — log(2;) {z is the parameter in the ice type model) passes from imaginary 
to real q^^'^ — + e~^. This shows a singularity in the coordinates but not 
necessarily the model itself at g = 4, n = ^/q = 2, so it seems that a role for 
surface roughening has not been excluded. 

If ric can be promoted to ~ 2.62 by surface roughening, then a 

Ge_3 = Go,2,/{p^ = DE?, might be available as the ground state space of an 
honest perturbation of i?o,3- Alternatively, the stability may at n = 2.62 be 
so slight as to be physically irrelevant. In either case, roughened or not, the 
quantum medium must certainly have topological dimension = 2 to admit 
anyons but the Hausdorff dimension of Y might approach 3. 

To be susceptible to the imposition of a topological symmetry (to force the 
system to be perpendicular to the ideal (p^+i)) we need Go,£ to be unstable 
(or in the f = 3 cases, nearly so). What will be the properties of a G^/ = 
Go/r\{p£+i)* if in fact such a ground state space is achieved? In § 2 the global 
properties of DE3 as a UTMF were discussed. These have implications for 
the local properties of a unit vector ^ e G^/ — Go,e n (pe+i)* and these 
contrast with a vector $ e Go,e. 

First, because ^ has the structure of a UTMF no information about 
the state ^ can be determined from observables acting on a disk D imbedded 
in Y, D G Y. Usual correlations such as (8) must be zero (or at least 
exponentially decaying) or else measurements in a large disk D (Z Y and 
"extrapolation" would reveal information about the state ^ on y . One might 
think with rapid decay of correlations, that observing \& in the classical basis 
we would see a "bubble phase" with a few tight global walls 7o C 7 forced by 
topology. But this is impossible such global lines could be locally detected 
on a disk D cY and would reveal information on ^ which is forbidden. In 
other words, a local operator could split the ground state whereas no local 
operator should effect more than an exponentially small splitting of G^^^. 

How is the paradox resolved? The domain walls "loops" in a typical 



48 



(according to L^— norm) component of \1/ will be very long probably space 
filling but at the same time not locally correlated. This behavior is seen 
already in the typical classical (i.e. observed) states of any toric code word 
[Kl], so the phenomena is not a surprise. 

Next we show that on (Y, A) the ground state space Go^i of ifo,^ is polylog 
extensively degenerate. Understanding this scaling is an important ingredi- 
ent in the perturbation theory. Fix a closed surface Y of genus {Y) = g and 
use the number of vertices v{A) as a measure of the combinatorial complex- 
ity of the triangulation A. Assume, for studying the v{Ai) := Vi — > oo 
asymptotics of Go, that the triangles of the triangulations Aj have bounded 
similarity type. This means that triangle shapes should not be arbitrarily 
distorted. In this regard "barycentric subdivision" is "bad" but more regular 
subdivisions are "good". 



Proposition 3.8. If genus (Y) = g > 1 then dim{Go,i{Y, A)) and dim(G';_^(F, A)) 
are O (f (A)^^~^). If genus (Y) = 1, the dimensions are 0{v{A)). 

Proof A pant is a 3— punctured 2— sphere. Fix a hyperbolic metric on 
Y. In a hyperbolic metric, d— isotopy classes have unique geodesic repre- 
sentatives. Using the Fenchel-Nielsen coordinates [Th] on a geodesic pants 
decomposition of (Y, A), we find, for each pant, 3 multiplicity parameters 
and 3 twist parameters each assuming 0{v^^'^) values, a definite fraction of 
which are mutually consistent. Since the Euler characteristic x (pant)= —1, 
there are 2g — 2 = |x(Y)| pants in the decomposition. The consistent pa- 
rameter settings define geodesic patterns, O {v^^~^) geodesic 1— manifolds 
compatible with A. These 1— manifolds are unique up to isotopy and have 
no trivial circles and so are also unique up to isotopy. But Go/ is defined 
as the perpendicular to the d— isotopy relation on 7i. □ 




bad 



good' 



Figure 3.2 
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Theorem 3.9. Assume A of Y is sufficiently fine. If {Oi} is a family of 
strongly local Hermitian operators on 7i with jgs X <z7i then there are three 
possibilities for X n Go/: 

(1) xnG'o,, = {0}, 

(2) X n Go/ = Go/, and 

(3) X n Go/ = Go/ n (p^+i)* = DEliY). 
The choice Oi — p\j^i realizes possibility (3). 

Note 3.10. By the Verlinda formulas dim DEi{Y), Y a closed surface, is 
asymptotically ^^^= (sin (the fraction of error converges to zero). 

This maybe compared to the much larger dim of Go/ (proposition 3.8). 

Proof: For each Oi there are (orthonormal) vectors fij spanning the k ten- 
sor factors on which Oi acts nontrivially so that Eo^i the lowest eigenspace of 
Oi, has the form span(/j i, /j 2, • • • /i.^J® (the remaining n — k factors) where 
5i — dim(£'o,i). The fij are assumed to be chosen coherently with respect 

to the natural isomorphism Oi = Cj. For each i © \fi,k){fi,k\ constitute 

k=l 

a skein relations Sj as explained above. Thus X fl Go/ consists of vectors 
orthogonal to the equivalence relation spanned by both d— isotopy and the 
skein relations {sj}. 

The easiest example is for the level= 1 theory where A = ie'^^^^, d = 
—A"^ — A^"^ = 1. In this case the Jones Wenzl projector p2 reads: P2 = )( — X 
or in combinatorial model: 




Figure 3.3 



Any combinatorial version of the smooth relation will give the same quo- 
tient and same Go/ H X provided the triangulation A is sufficiently fine. 
When d = 1 generalized isotopy simply means isotopy and deletion of circles 
bounding disks; adding Si = yields unoriented Z2— homology — as shown 
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in Figure 3.3 — as the quotient equivalence relation. So, for example on a 
closed surface Y, the space X fl Go/ has dimension 2^^*^^^ bi being the first 
Betti number, and is identified with functions Hi[Y; Z2) — > C. 

For any skein relation s the subset < Si > (iGo/ C TC, each Si specializing 
s at locations Ui, is an ideal of the surface algebra Go.e as defined in section 
2. The uniqueness Thm 2.5 implies < Si >= Go/, {0}, or {pe+i) (the latter 
occurs when Si represents pe+i on Ui, i.e, Si = p^+i). These correspond 
respectively to the three alternatives in the theorem. □ 

Observation 3.11. Although eV = e T^a^, is promising perturbation to find 

2^ 

DEi, its ground state vector 9o = 2^2^^ (— 1)*'~^^A: has exponentially large 
Ho/— expectation, i > 1, 

{0o\Ho/\eo) >e^'^ (3.15) 
some a > 0, and L the refinement scale. 

Proof: If a spin configuration is chosen uniformly at random from the 2" 
possibilities it follows from an easy independence argument that the mean 
number of circles # of that configuration is 0{v) and the standard deviation 
is s.d. (#) > (9(f ^/^), V the number of vertices of A(y). For d > 1, line 3.15 
is deduced using the "circle" term of Ho/ together with the above inequality 
on standard deviation. □ 
To consider the possibility of frustration, which is outside our algebraic 
ansatz (but quite possible in a fundamental description of an anyonic system, 
see § 4, (2) Uniqueness) we should think about how a general local operator 
could act on Go/. Such (non scalar) actions are possible precisely because Go/ 
is not a topological modular functoiQ (there is no disk axiom here) so a local 
operators such as O may detect statistical information on the topology of a 
state $ G Go/- For example the presence of a bond in the domain wall 7 may 
have polynomial influence on a global topological event. This phenomena is 
familiar from percolation: The state of a single bond in the middle of an 
i X i piece of lattice has influence, at Pc, on the existence of a percolating 
cluster joining two opposite boundary segments which decays as (for the 
triangular lattice) [LSW] and [SW]. O may determine an effective reduction 
-Eao C Go/ based on some local statistical difference between the random 

''For £ < 2 Go/ presumably has gapless modes which are un - topological. For £ > 2 
string tension allows local measurements to yield a global inference. 
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components of different G 7i which the local operators |si)(si| detect. 
The sites of these operators should be located and oriented randomly w.r.t. 
the domain wall 7^ of $fc so it is reasonable to believed that the deviation 
of O from scalar x identity will be attributable to the topologically essential 
part 7^ of 7fe since other features may be lost in averaging over the sites. For 
example, 7^ may appear "straighter" than the "foam" 7^X7^. That is, O 
may favor or disfavor a topologically complex essential-domain-wall 7]^ over, 
say, a simpler essential-domain- wall 7^. 

If complexity is favored the ansatz is still capable of describing the fi- 
nal reduced ground state space, though for an indirect reason: A reduction 
Ex^ C Go,e, or a series of reductions, which preserves the complex topo- 
logical representatives will retain many representatives for each class in the 
quotient Go// (pe+i) modular functor so topological information will not be 
lost, Ex^ n ipi+i)* = Go/H (pi+i)*. (But note that the intermediate subspace 
Ex^ fails to respect the multiplicative/tensor structure on Go/.). 

If topological complexity is disfavored in Ex^ , then the action of O could 
destroy the topological information of the modular functor by killing all 
classes except for the simplest, 7"'" = 0, corresponding to no essential domain 
wall and foam covering all oi Y. In this case the ansatz is not applicable. 
But fortunately the sign of e in the perturbation eV may be possible to con- 
trol, placing the system (7Y, A, H^) into the favorable regime, where a further 
reduction, perhaps at, higher order, could still find the modular functor. 

Because of the topological character of skein relation, the final quotient 
Go,i/{p£+i) is the same regardless of which sites and with what coefficient 
norm the various pl_^_i relation is enforced. Applying p^+i requires no ho- 
mogeniety of input and in fact "smooths" local disturbances. 

4 The evidence for a Chern-Simons phase 

There is no proof or "derivation " of a stable (gapped) phase 
G^/ = Go// < Pi+i >= DEi but there is evidence in the form of anal- 
ogy and internal consistencies under the headings: "UTMF" , "uniqueness" , 
and "positivity" . The first two have been discussed and are only now sum- 
marized, the third is presented in more detail. 

1) UTMF. The quotient algebra Go// < Pe+i >= DEi is an anyonic 
system, in mathematical terms a UTMF. This in itself is consistency check: 
The presence of an anyonic quotient system. The spatial correlation scale for 
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topological information is zero so we expect a gap protecting the topological 
degrees of freedom. The system is the quantum double of a Chern-Simons 
theory and experience with the FQHE as a Chern-Simons theory has prepared 
us to believe that these beautiful structures can self-organize in nature from 
the simplest underlying Hamiltonians - e.g. Coulomb repulsion in a Landau 
level. 

Mathematically the double has the form of the algebra of operators on 
some (fictitious) FQHE - like system "X" : domain walls realize the Wilson 
loop operators on X. The double is freed from chiral asymmetry and the 
extreme physical conditions required to break time reversal symmetry. The 
double is a better place to look for a reahzation protected by a large spectral 
gap. 

2) Uniqueness. There is a unique candidate model G^^/ < Pi+i > re- 
specting the local or "multiplicative" structure of Go,€ (Thm3.9). Uniqueness 
suggests that there is a sharp boundary: there are no slightly larger quotients 
which could include low frequency excitations. Thus the simplest expecta- 
tion, that the reduction G^^£,i > 2, be one dimensional (non-degenerate), 
cannot be achieved as a joint ground state jgs of local projectors but requires 
frustration . A subtle point is involved here. The most interesting candidates 
in solid state physics for topological states are highly frustrated systems 
when written out in their fundamental degrees of freedom. This does not 
mean that they cannot have effective descriptions as a jgs (= unfrustrated) . 
In fact a very general topological argument suggest that a phase with the 
structure of a TQFT always will. In local models (e.g. [TV]), product struc- 
tures Y X I yield projectors H{Y) — ^ whose image is the underlying 
UTMF, V{Y). However by building Y x I from overlapping local product 
structures, p can be factored into a commuting family of local projectors {pi} 
so V(Y) = image p = jgs{pi}. 

Finally, uniqueness creates an aesthetic bias: Could nature really turn 
down such a possibility? 

We now turn to the final consistency check. 

3) Positivity of the Meirkov trace pairing on Go,e/{pe+i)- 

On a surface Y with or without boundary, the Hilbert space Go//{pe+i) = 
Go/ n (pe^i)* inherits a Hermitian inner product (, )geom. by inclusion in 
7i, the space of all spin configurations. (Beginning with {|+), | — )} as an 
orthogonal basis for C^, the Hilbert space H acquires the tensor product 
pairing which may be restricted to fl (p^+i)*.) On the other hand. 
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G'^^/{pi+i) = DEi{Y) has a topologically defined "Markov trace" Hcrmi- 
tian inner product ( , )top. corresponding to its structure as a UTMF (see 
Definition 4.2.). 

Proposition 4.1. Suppose that the perturbed Hamiltonian H^ i has a spec- 
tral gap above its ground state Go/ H (pe+i)* , then up to a correction which 
is exponentially small in the refinement scale of the triangulation A, ( , )geom. 
and ( , )top. arc proportional: ( , )gcom. — c( , )top. for some real number c 7^ 0. 

The definition of ( , )top. is recalled below. 

Definition 4.2. Extending the definition given in section 2: If 71 and 72 

are domain walls in Y with identical boundary data then 7 = 71 U 72 defines 
a link inYxS^ where Y is Y with its boundary capped by disks. (Place 71 
and 72 on disjoint 6'— levels 9i and 62, then bend paired endpoints to meet at 
the intermediate level (^^i + 62) /2 w.r.t. the orientation.) Regarding 7 as 
labelled by the 2— dimensional representation, (71,72) := Witten invariant 
(y X S\'y) at level £, see [Wi], [RT], and [BHMV]. We refer to this pairing, 
also, as the Markov trace pairing. 

The combinatorial properties the code space C := Go/ n (pe+i)* are such 
that local operators (in fact any operator supported on some topological disk 
D G Y) cannot extract or modify information in C. (It is true that a local 
operator can rotate C to C, C _L C, but nondestructive (of details within 
C) measurements allow the error to be corrected by a physical operator 

acting on 7Y so that the composition C — > C — > C is the identity idc-) 
The code property is a kind of combinatorial/topological rigidity and it is 
quite natural that, if achieved in a ground state space, that space should be 
protected by a spectral gap. 

As we argue for proposition 4.1, a final consistency check emerges: 

gap + code =^ positivity of Markov trace pairing. (4.1) 

This explains how the Markov pairing is picked out and why its (indefi- 
nite) Galois conjugates are unreahzable as stable phases. The Markov trace 
pairing is known to be positive [J], [FN WW] precisely for our choice of ^4, 

A = ie'^*/^'', d = —A? — A~'^, and it being topologically defined is automat- 
ically invariant under the mapping-class-group of {Y,dY). For other roots 
of unity A the resulting Hermitian pairings are of mixed sign so cannot, by 
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positivity of ( , )geom. on Ti and proposition 4.1, correspond to stable physical 
phases. For example, at level £ = 3, Galois conjugate choice d' = is not 
a physical. No ground state space modelling / < p^' > could have a gap. 
The corresponding "UTMFs" are only "unitary" with respect to mixed (p, q) 
Lorentz form [BHMV], constructed for each labelled surface of the theory. 
These structures cannot be induced by restricting the standard Hermitian 
pairing ( , )geom. on U. 

By choosing d = 2 cos ^ as we have ensured tr (a, a) is positive. If G^/ 
has a gap and is naturally identified with Go// (p^+i) then 4.1 implies posi- 
tively of the Markov trace. Positivity is demonstrated by showing that, up to 
an error exponentially small in the refinement scale L, that the Markov trace 
is in the similarity class of the Hermitian form induced from the standard 
inner product on 7i. This is the argument; it is not mathematically rigorous 
as the "imaginary - time = space ansatz" is employed, but we hope that is 
convincing physically. 

Argument, Proposition 4.1. A surface (F, A) can be gradually changed 
by bringing bonds in and out of the triangulation (and perhaps adding or 
deleting vertices). With patience, a Dehn twist can be effected. This takes 
Oir?^ moves on an n x n square grid torus T^. Similarly a braid generator 
for quasiparticle excitations on a disk takes 0{r?^ such moves where n is 
the number of bonds in a loop surrounding the two quasiparticles. These 
changes can eventually return A to a homeomorphic, though now twisted, 
image of itself. 




before after 
Figure 4.1 



If ifg has a gap, bounded as we change A (on which depends), the 
adiabatic theorem will define, in the slow deformation limit: deformation 
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speed << gap, a time evolution of vectors in ^ C Ti, ie[0, 1]. At time 
t = 0,Gt:^o = and finally at time t = 1, G^ i = G^ again. This evolution 
is (incidentally) identical with the one induced by the canonical connection 
on the universal topological bundle of {/c - plane, vector in k - plane } — > 
{k - planes}. Prom the assumption of a gap= 5, one can argue that this 
monodromy for a Dehn or braid twist is accomplished by a composition of 
0{6~^n'^) local operators, or more precisely operators At which have only 
an exponentially small nonlocal part. This means that for Pauli matrices 
at sites i and j, the commutator satisfies: || [al.At,a^] ||< Co e"'^^"*"-^" for 
all indices x,y,i,j,t, and some positive constants Co and ci. We call such 
operators quasi-local. The essential point is that the local disturbance caused 
by modifying Ht near a bond to Ht+i dies away exponentially in imaginary 
time and hence in space. Let us ignore the exponential tail (which will lead to 
a manageable error term), and think of the monodromy as a composition of 

0(5"^n^) local operators: monodromy = H^t"'^- ^^^h i in a discretized 

t 

unit interval with 0{d~^n^) points, Al°^ is a local unitary operator G Hom 
(7i, H) which carries the code subspace Ge,t of H at time — t to the code 
subspace Ge,t+i at time — t + 1. 

Adiabatic evolution has provided us with one (local) representation, Pgeom. '■ 
TTi (moduli space) — > -Pf^gcom.(G'e.t), from the fundamental group of moduli 
space (F) (in our discrete context moduli space is the space of triangulations 
of Y) to the projective unitary transformations of the perturbed ground 
state space. The subscript geom. signifies that PU is defined with respect 
to (, )geom.- On the other hand, assuming a spectral gap above G^^t- there 
is a physical argument that a second, topologically defined, representation 
is also local. This representation: ptop. '■ tti (moduli space) — > -Pt^top. (G^e,t) 
is defined into the projective unitaries w.r.t. ( , )top. by deforming the tri- 
angulation while leaving the formal picture (:= superposition of domain 
walls) topologically invariant. This representation can be defined by choos- 
ing a local rotation which interpolates between the conditions (that define 
Go,t n ipe+i)*) in force at time = t but not t+1 and those in force at time 
— t + 1 but not t. What is not immediate is whether the effect of this lo- 
cal rotation on the jgs can be achieved by an operator Af on Ti. which is 
quasi-local. But the existence of a quasi-local At can be argued based on 
the "imaginary — time = space" ansatz (§ 3 lines (8)-(14)). Similarly, if we 
view the ground state Ge,t as a local excitation of i?e,t+i, but one without 
topological content, we expect that they can be annihilated by a quasi-local 
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A- 

But if a local operator carries one code space into another, that oper- 
ator restricted to the first code space is unique, up to a scalar, among all 
restrictions of such local operators. This is particularly clear in the present 
case when the operators are unitary and all the code spaces have the same 
dimension. Suppose both A and B are unitary operators carrying Ci into 
C2, then B'^ o A\c-^ : Ci — > Ci is also local and so multiplication by some 
unit norm scalar A. Thus B\c^ = 

So assuming a gap, the proceeding observation shows first that Pgeom. and 
Ptop. are both actually well defined as maps from the fundamental group (see 
Figure 4.2) and second that Pgeom. will be projectively the same ptop. up to 
an error exponentially small in the refinement scale L (when measured in 
the operator norm). The latter, ptop is simply parallel transport in Witten's 
[Wi] projectively flat connection on the modular functor bundle V{Yt) over 
the moduli space of surfaces {Yt} {t now an arbitrary parameter). Projec- 
tive flatness as well as uniqueness of this connection follow formally from 
locality properties: As the surface is gradually changed (discretely this is 
done by moves on the triangulation A) the two surfaces Yt and Yt+i can be 
canonically identified in the complement of a disk D supporting the changing 
bonds, and the identification can be extended arbitrarily over D. Prom the 
disk axiom and the gluing axiom of section 2, we have a unique canonical 
projective isomorphism of modular functors V{Yt) — > ViYt^i). This deter- 
mines, via differentiation, a unique connection. Projective flatness follows 
by applying this uniqueness to a loop of identifications representing a small 
cycle of changes to A collectively supported in a disk D <zY. Similar loops 
span the relations in tti (moduli spaces). 
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Above, we see atypical 5 -cycle relation induced by sweeping 
out the boundary of a 4-simplex with a triangle. 



Figure 4.2 

Let V = DEi £ 7^ 1, 2, or 4. It is known [FLW2] that for a sphere with 
4 or more punctures (or a higher genus surface) the braid (or mapping class 
group) acts densely in the projective unitary transformations of each label 
sector PUtop.{V{Y, t )). But identifying V{Y, t ) with a ground state space 
G^/ (followed by the idempotent t defined in §3), we have on the one hand 
the adiabatic evolution which must be unitary w.r.t. the Hermitian pairing 
induced from the standard Hermitian pairing on Ti, and on the other hand, 
exponentially close to this, transport in Witten's connection. Both define 
(nearly) the same dense homomorphism from tti := the fundamental group 
of moduli space: 

Pgeom. Ptop. : VTi End{V{Y, t)) . (4.2) 

In the case Y is planar and all boundary labels equal, tti is a familiar braid 
group. 

It follows from the rapid approximation algorithm [KSV],[So],[K2] of ele- 
ments of PUgeom., PUtop. C End(y{Y, t )) by words in tti, that the induced 
Hermitian metric on V must be exponentially close (in L) to the intrinsic 
UTMF metric on V up to the overall scalar c. The mathematical fact that 
we are using here is that Hermitian pairings can be recovered from their 
symmetries: 

Lemma 4.3. Suppose a vector space V has Hermitian (but not necessar- 
ily positive definite) inner products (, )i and (, )2 with symmetry groups 
U11U2 C End(\^) respectively, if Ui — U2 then { , )i — c{ , )2 for some real 
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constant c ^ 0. Furthermore if Ui ^ Ui but instead if for all Ai G Ui there 
exists an A2 G U2 with \ \Ai — < e > 0, and for all A^ G U2 there exists 
an Ai G Ui with ||yl2 — ^i|| < e, then for all v G V, {v,v)l/ {v,v)l = const. 

Proof: Up to linear conjugacy the type of the form is determined by dimen- 
sion, signature and nuUity. If Ui = U2 (even approximately) these invariants 
agree. Let MeEnd(l^) transform and {,)2 

{{M^v,Mw)i = {v,w)2). Then U^M = MU2 so if Ui = U2 =: U then 
M normalizes U. Let s + t = dim{V). In PGL{s,t;C), PU{s,t) is its 
own normalizer establishing the lemma when U is nonsingular. If the forms 
have radicals, these must agree and the preceding argument applies mod- 
ulo the radical. Finally, in the case Ui and U2 are not identical but have 
Hausdorff distance = e, a counter example to the lemma would yield a Lie 
algebra element a G pgi{s,t)\j)u(s,t) with ada{pu(s,t)) C pu{s,t), but 
pu{s,t) C pgi{s,t) is a maximal proper sub Lie algebra. □ 
Consequently, G^^i is not just linearly V but metrically V, provided H^ ^ 
has a gap. 

The spectral gap assumption implies that the combinatorically defined 
Markov trace pairing is induced from the standard inner product of Ti. The 
Markov trace pairing is a rather intricate structure in its relation to gluing 
(see axiom 2). That it arises from a simple assumption can be viewed as a 
valuable "consistency check" on that assumption - the existence of a spectral 
gap above Ge,e. 

5 The i^^e/i medium as a quantum computer 

In the literature one finds at least three polynomially equivalent models of 
quantum computation defined: g— Turing machine [D], g— circuit model [Y], 
5— cellular automata [LI]. Nearly all proposed architectures ([NC] is an ex- 
cellent survey) presume localization of the fundamental degrees of freedom. 
This may be called the "qubit approach" although qunit might be more pre- 
cise since there is nothing special about two state systems, the number n of of 
states per site may even be infinite, as in optical cavity models — what is im- 
portant in these architectures is the tensorial structure of the computational 
degrees of freedom. 

However, there is another approach [FKLW] in which the global tensor 
structure becomes redundant. The physical degrees of freedom still have a 



59 



local tensor structure — as is universal in quantum mechanics ^ but these 
are never touched directly. Instead a system is engineered so that these 
local degrees interact through a Hamiltonian H whose eigenstates E\ are 
highly degenerate codes spaces capable of storing, protecting and processing 
quantum information. For us Ex will be the internal symmetries of an anyonic 
system in which position coordinates have been frozen out. The processing 
will consist of braiding anyons in (2 + 1)— dimensional space time. 

To make sense of this, consider a definition of " universal quantum com- 
puter" which does not presuppose any tensor decomposition. We need: 

1. h: A Hilbert h space on which to act. (Its dimension should scale 
exponentially in a physical parameter.) 

2. We need to be able to initiahze the system. 

3. p: Operations — > U{h), a representation of some group (or at least 
semigroup) of operation on the unitary transformations of h which can 
be physically implemented - preferably with error scaling like e~'^°'^^*^'^* ^ 
for some physical parameter L. (Lack of such scaling is the Achilles 
heel of qubit models.) The representation p should have dense image 
in SU{h): This, together with the rapid convergence property of dense 
subgroups of U (h) , ensures universality. 

4. Compiler: This is a classical computer which takes a g— algorithm and 
an instance, e.g. Shor's poly time factoring algorithm [S2] and a thou- 
sand bit integer, and maps the pair into a string s of operations as in 
(3.). 

5. ^f-. The result of the quantum portion of the calculation is a final state 

= p(s)*o. 

6. Observation: There must be a Hermitian operator which serves as the 
observation: projecting into an eignstate with probability \ax\'^, 
^f — Eoa^a- The eigenvalue A is what is actually observed. 

7. Answer: Another poly-time classical computation is now made to con- 
vert the observed eigenvalue, perhaps for many executions of 1) — > 6), 
into a probabilistic output. The class of problems that can be answered 
in polynomial time by 1) — > 7) with bounded error probability (say 
error < |) is called BQP. For example, factoring [S2] is in BQP. 
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Computer scientists believe, and cryptographers hope, that factoring 
is not in the corresponding classical computational class BPP. 

The reader can easily take any qubit architecture (see [NC] for details 
of these) and fit it into the proceeding format. Let us now do this for our 
anyonic system with Hamiltonian H^ £. As explained in the introduction, the 
system was chosen to have two spatial dimensions topologically, i.e. to live on 
a triangulated surface (F, A), so that exotic statistics become a possibility. 
By mathematical excising neighborhoods of the excitations we reduce to the 
case of the studying the ground state space Gt of a time dependent H^^i^t, 
on a highly punctured surface with labelled boundary. The subscript t 
reminds us of the time dependence of the surface A) as the position 
of the punctures evolves. The ground state space Gf describes the internal 
symmetries of a collection of quasiparticle (anyon) excitations whose spatial 
locations are t dependent. The space Gt„ is a representation space for the 
braid group (or generalized braid group) which describes the motion of these 
quasiparticles on Y. Because of the presumed spectral gap, the quasipar- 
ticles are expected to have exponentially decaying tails. Chopping off and 
ignoring these tails amounts to the puncturing Y at the quasiparticles. This 
identifies the excited state j e Ex^t containing anyons on Y with a ground 
state \E'o,t € Gt on a multiply punctured Y~ with labeled boundary. So by 
puncturing and labeling the surface Y, a ground state in Gt can be used to 
represent the anyonic state e E^^t, so the discussion of § 3 applies to Ex^f 

Let us walk through steps 1 through 7 for our anyonic model, though it 
is not efficient to do this in strict order. 

1. & 2. /i = Ex,t = Gt: In [FLWl] and [FKLW] abstract anyonic models for 
(but with no known Hamiltonian H) were analyzed algebraically. In 
[F2] an explicit but artificial Hamiltonian was given as an existence 
theorem. The UTMFs of [FLW1],[F2] and [FKLW] required a two 
(with care 1.5) quasiparticle pairs per qubit simulated. DEi is a closely 
related UTMF and for £ = 3 a similar encryption yields one qubit 
per 1.5 pairs of (0,2) type excitatitons. Physically one imagines a 
disk of quantum media governed by Ht and trivial outer boundary 
condition, lying in its (nondegenerate) ground state = "the vacuum" . 
The steps required to build h are a subset of those discussed in [BK] 
in connection with their CS2 model. The disk is struck in some way 
(with a hammer?) at a point to create a pair of excitations. Already 
in building Ex,t we need measurement to tell if the newly created pair 
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is type ((0,2); (0,2)). If the pair is of this type, we keep it, if not it 
is returned to the vacuum. Repeat (perhaps thousands of times) until 
a sufficiently large Hilbert space Ex^t{Y) — Gt(Y^) = h, and initial 
vector, ^'o G Gt is realized. The initial state ^'o is determined by 
the condition that all circles surrounding (not separating) the created 
pairs acquire label (0,0). How many pairs are required depends on 
the problem instance. For example, for factoring problem, it is a small 
multiple of the number of bits of the number to be factored. 

6. Mecisurement: The creation process is probabilistic. So even at the 

start, there must be a local observation which tells us which anyon pairs 
have been created ((0, 0), (0, 0)); ((0, 2), (0, 2)); (2, 0), (2, 0); or ((2, 2), (2, 2)). 
One hopes that will not require a "topological microscope". Because 
quasi particles are arrangements of elementary degrees of freedom spins, 
charges, etc. . . of the system, one expects each quasiparticle when 
examined electromagnetically to have its own unique signature: e.g. 
quadruple moment etc. ... In this view, localized quasiparticles would 
always be "measured" by their environment and never lie is superpo- 
sitions. However, it is essential for quantum computation that a well 
separated pair of quasiparticles — before being fused — could be in a 
super position of collective states. Another idea [SF], discussed in the 
introduction, is that a phase transition be employed for measurement. 

3. Braiding: The group of operations is the braid group of quasiparticles 
moving on the disk. In order to implement this mathematically known 
representation on h we need to be able to grab hold of the quasiparticle 
and, within some allowable dispersion corridor, nudge it along to exe- 
cute the braid s dictated by the output of the classical compiler, step 4. 
This, like observation, should in principle be possible, using the charac- 
teristic electric and magnetic attributes of the nontrivial quasiparticles 
(whatever they are eventually measured to be). It may be possible to 
design wells that trap, and when desired, move specific quasiparticles. 

5. ^/ is the internal state after braiding ^/ = p(s)^o- 

6. Observation: We already discussed the necessity to observe halves of 
newly created pairs. To read out quantum information after braiding 
take two quasiparticles in the system ^ / and fuse them. Although they 
will retain their individual identities during braiding and still be both 
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of type (0,2), after fusing, two outcomes are possible: (0,0) or (0,2). 
The probabilities attached to these outcomes is the classical distillation 
of quantum information equivalent to measuring a qubit in the usual 
architecture (see [FKLW] for details of the read out and its relation 
to quantum topology and the Jones polynomial.) The braiding has 
rearranged, in an exponentially intricate fashion the structure of the 
composite pairs of (0,2)— quasiparticles. This recoupling is the heart 
of the computation. 

7. &c 4. The final conversion of eigenvalues observed to a probabilistic output is 
the same as for the qubit architecture. The structure of the compiler is 
also similar but must include a rapid approximation algorithm [K2] [So] 
subroutine. 

We have presented -f^e,3,t as a theoretical candidate for an anyonic medium 
capable of universal quantum computation. Its experimental realization 
would a landmark. 

6 Appendix. Ideals in Temperley-Lieb Cater- 
gory 

Frederick M. Goodman and Hans Wenzl 

This appendix contains a proof of the following result, which is used in the 
paper of Michael Freedman, A magnetic model with a possible Chern-Simons 
phase. 

Theorem 0.1. When the parameter d is equal to 2 cos(j7r/n) with n > 3 and 
j coprime to n, then the Temperley-Lieb category has exactly one non-zero, 
proper ideal, namely the ideal of negligible morphisms. For all other values 
of d, the Temperley-Lieb category has no non-zero, proper tensor ideal. 

We are grateful to Michael Freedman for bringing the question of tensor 
ideals in the Temperley-Lieb category to our attention and for allowing us 
to present the proof as an appendix to his paper. 

Our notation in the appendix differs slightly from that in the main text. 
We write t instead of —A^, Tn for the Temperley-Lieb algebra with n strands, 
and TL for the Temperley-Lieb category. We trust that this notational vari- 
ance will not cause the reader any difficulty. 

This appendix can be read independently of the main text. 
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1 The Temperley-Lieb Category 



1.1 The Generic Temperley Lieb Category 

Let t be an indeterminant over C, and let d = The generic Temperley 

Lieb category TL is a strict tensor categor whose objects are elements of 
No = (0, 1,2,...}. The set of morphisms Hom(m, n) from m to n is a C{t) 
vector space described as follows: 

If n — m is odd, then Hom(m, n) is the zero vector space. 

For n — m even, we first define (m, n)-TL diagrams, consisting of: 

1. a closed rectangle R in the plane with two opposite edges designated 
as top and bottom. 

2. m marked points (vertices) on the top edge and n marked points on 
the bottom edges. 

3. (n + m)/2 smooth curves (or "strands") in R such that for each curve 
7, (97 = 7 n dR consists of two of the n + m marked points, and such 
that the curves are pairwise non-intersecting. 



Two such diagrams are equivalent if they induce the same pairing of the 
n + m marked points. Hom(m, n) is defined to be the C{t) vector space with 
basis the set of equivalence classes of (m, n)-TL diagrams; we will refer to 
equivalence classes of diagrams simply as diagrams. 




Figure 1.3: A (5, 7) -Temperley Lieb Diagram 
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The composition of morpliisms is defined first on the level of diagrams. 
The composition ba of an (m, n)-diagram b and an {i, m)-diagram a is defined 
by the following steps: 

1. Juxtapose the rectangles of a and b, identifying the bottom edge of a 
(with its m marked points) with the top edge of b (with its m marked 
points) . 

2. Remove from the resulting rectangle any closed loops in its interior. 
The result is a (n, ^)-diagram c. 

3. The product ba is d^c, where r is the number of closed loops removed. 

The composition product evidently respects equivalence of diagrams, and 
extends uniquely to a bilinear product 

Hom(m, n) x Hom(^, m) — > Hom(^, n), 

hence to a linear map 

Hom(m, n) (8) Hom(£, m) — > Hom(£, n). 

The tensor product of objects in TL is given hj n ® n' = n + n' . The 
tensor product of morphisms is defined by horizontal juxtposition. More 
exactly, the tensor a 6 of an (n, m)-TL diagram a and an {n', m')-diagram 
b is defined by horizontal juxtposition of the diagrams, the result being an 
(n + n', m + m')-TL diagram. 

The tensor product extends uniquely to a bilinear product 

Hom(m, n) x Hom(m', n') — > Hom(m + m', n + n'), 

hence to a linear map 

Hom(m, n) <S> Hom(m', n') — > Hom(m + m', n + n'). 

For each n e No, T„ := End(n) is a C(t)-algebra, with the composition 
product. The identity 1„ of T(n) is the diagram with n vertical (non-crossing) 
strands. We have canonical embeddings of T„ into Tn+m given by x i— > x<S>lm- 
lim > n with m — n even, there also exist obvious embeddings of Hom(n, m) 
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and Hom(m, n) into T,„ as follows: If fl and U denote the morphisms in 
Hom(0, 2) and Hom(2, 0), then we have hnear embeddings 

a e Hom(n, m)^a® u®^"*-")/^ e 

and 

h e Hom(m, n)^h® e T^. 

Note that these maps have left inverses which are given by premultiph- 
cation by an element of Hom(n, m) in the first case, and postmultiplication 
by an element of Hom(m, n) in the second. Namely, 

a = d-("*-")/2(a (g) U®("^-")/2) o (1„ ® f^^{m-n)/2^^ 

and 

b = d~^"'~''^/^{ln u^i"'-^)/^) o (6 (g) n®("^-")/2) 

By an ideal J in TL we shall mean a vector subspace of 0^ ^ Hom(n, m) 
which is closed under composition and tensor product with arbitrary mor- 
phisms. That is, if a, b are composible morphisms, and one of them is in J, 
then the composition ah is in J; and if a, b are any morphisms, and one of 
them is in J, then the tensor product a (g) 6 is in J. 

Note that any ideal is closed under the embeddings described just above, 
and under their left inverses. 

1.2 Specializations and evaluable morphisms. 

For any r e C, we define the specialization TL(r) of the Temperley Lieb 
category at r, which is obtained by replacing the indeterminant thy t. More 
exactly, the objects of TL(r) are again elements of Nq, the set of morphisms 
Hom(m, n) (r) is the C-vector space with basis the set of (m, n)-TL diagrams, 
and the composition rule is as before, except that d is replaced by (i(r) = 
(r + T~^). Tensor products are defined as before. T„(r) :— End(n) is a 
complex algebra, and x i— > x®!^ defines a canonical embedding of ^^(t) into 
T„+^(t). One also has embeddings Hom(m,n) — > T„ and Hom(n, m) — > T^, 
when m < n, as before. An ideal in TL(r) again means a subspace of 
m Hom(n, m) which is closed under composition and tensor product with 
arbitrary morphisms. 

Let C{t)r be the ring of rational functions without pole at r. The set of 
evaluable morphisms in Hom(m, n) is the C(i)T--span of the basis of (n, m)- 
TL diagrams. Note that the composition and tensor product of evaluable 
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morphisms are evaluable. We have an evaluation map from the set of evalu- 
able morphisms to morphisms of TL(t) defined by 

where the Sj are in C(i)^, and the aj are TL-diagrams. We write x 
x{t) for the evaluation map. The evaluation map is a homomorphism for 
the composition and tensor products. In particular, one has a C-algebra 
homomorphism from tlie algebra of evaluable cndomorpliisms of n to the 
algebra T„(t) of endomorphisms of n in TL(t). 

The principle of constancy of dimension is an important tool for analyzing 
the specialized categories TL(t). We state it in the form which we need here: 

Proposition 1.1. Let e & Tn and f e be evaluable idempotents in the 
generic Temperley Lieb category. Let A be the C(t)-span in Hom(m, n) 
of a certain set of {m,n)-TL diagrams, and let A{r) be the C-span in 
Hom(m, n) (r) of the same set of diagrams. Then 

dimc(t) eAf = dime e(T)A(T)/(T). 

Proof. Let X denote the set of TL diagrams spanning A. Clearly 

dimc(t) A = dime A(r) = \X\. 

Choose a basis of e(r)A(r)/(r) of the form {e(r)a;/(r) : x G ^o}) where 
Xq is some subset oi X. If the set {exf : x e Xq}, were hnearly dependent 
over C{t), then it would be linearly dependent over C[t], and evaluating at 
T would give a linear dependence of {e{T)xf{T) : x e Xq} over C. It follows 
that 

dimc(t) eAf > dime e(r)A(r)/(r). 

But one has similar inequalities with e replaced by 1 — e and/or / replaced 
by 1 — /. If any of the inequalities were strict, then adding them would give 
dimc(t) A > dime A(t), a contradiction. □ 

1.3 The Markov trace 

The Markov trace Tr = Tr^^ is defined on T„ (or on T„(t)) by the following 
picture, which represents an element in End(O) = C(t) (resp. End(O) = C). 
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= Tr(a) e End(O) 



Figure 1.4: The categorical trace of an element a e T„. 

On an (n, n)-TL diagram a G T„, the trace is evaluated geometrically 
by closing up the diagram as in the figure, and counting the number c(a) of 
components (closed loops); then Tr(a) = d''^"'\ 

It will be useful to give the following inductive description of closing up 
a diagram. Wc define a map e„ : T„+i Tn (known as a conditional expec- 
tation in operator algebras) by only closing up the last strand; algebraically 
it can be defined by 

a e T„+i 1-^ (l„(8)U)o(a(8)l)o(l„(8)n). 

li k > n, the map en,k is defined by en,k = £n ° £n+i ■■■ ° Sfc-i- It follows 
from the definitions that Tr(a) = £o,n for ^ £ Tn. 

It is well-known that Tr is indeed a functional satisfying Tr(a6) = Tr(6a); 
one easily checks that this equality is even true if a G Hom(n, m) and b G 
Hom(m, n). We need the following well-known fact: 

Lemma 1.2. Let f G Tn+m (^nd let p G T„ such that (pCg> lm)/(p® Im) = /? 
where p is a minimal idempotent in T„. T/ien £n,n+m(/) = IP, where 7 = 
Trn+m{f)/Trn{p) 

Proof. It follows from the definitions that 

Pen,n+m{f)P = £n,n+m((P ^ lm)/(P 1^)) = £n,n+m(/)- 

As p is a minimal idempotent in Tn, Sn,n+m{f) — IP, for some scalar 7. 
Moreover, by our definition of trace, we have Trn+m{f) — Trn{Sn,n+m{f)) — 
jTrn{p)- This determines the value of 7. □ 
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The negligible morphisms Neg(n, m) are defined to be all elements a G 
Hom(n,m) for which Tr{ab) = for all b G Hom(m, n). It is well-known 
that the set of all negligible morphisms form an ideal in TL. 



2 The structure of the Temperley Lieb alge- 
bras 

2.1 The generic Temperley Lieb algebras 

Recall that a Young diagram A = [Ai,A2, .... A^] is a left justified array of 
boxes with Aj boxes in the i-th row and Aj > Aj+i for all i. For example, 



[5, 3] = m 



n 



All Young diagrams in this paper will have at most two rows. For A a Young 
diagram with n boxes, a Young tableau of shape A is a filling of A with the 
numbers 1 through n so that the numbers increase in each row and column. 
The number of Young tableax of shape A is denoted by fx- 

The generic Temperley Lieb algebras T„ are known ([2^) to decompose 
as direct sums of full matrix algebras over the field C(t), T„ = 0;)^Ta, where 
the sum is over all Young diagrams A with n boxes (and with no more than 
two rows), and Tx is isomorphic to an fx-^J-fx matrix algebra. 

When A and fi are Young diagrams of size n and n + 1, one has a (non- 
unital) homomorphism of Tx into given by x i— >■ {x ^ l)Zfj_, where 
denotes the minimal central idempotent in T„+i such that = Tn+iz^. Let 
gx,fj. denote the rank of {e ® l)z^, where e is any minimal idempotent in Tx- 
It is known that gx^^ = 1 in case /i is obtained from A by adding one box, 
and gx,f_i = otherwise. 

One can describe the embedding of T„ into T„+i by a Bratteli diagram 
(or induction-restriction diagram), which is a bipartite graph with vertices 
labelled by two-row Young diagrams of size n and n + 1 (corresponding 
to the simple components of T„ and T„+i) and with gx^^ edges joining the 
vertices labelled by A and /i. That is A and /i are joined by an edge precisely 
when /i is obtained from A by adding one box. The sequence of embeddings 
To Ti T2 — > . . . is described by a multilevel Bratteli diagram, as shown 



in Figure 2.5 



A tableau of shape A may be identified with an increasing sequence of 
Young diagrams beginning with the empty diagram and ending at A; namely 
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the j'-th diagram in the sequence is the subdiagram of A containing the num- 
bers 1, 2, . . . , j. Such a sequence may also be interpreted as a path on the 
Bratteh diagram of Figure 2.5, beginning at the empty diagram and ending 
at A. 

2.2 Path idempotents 

One can define a famihy of minimal idempotents pt in T„, labelled by paths 
t of length n on the Bratteli diagram (or equivalently, by Young tableaux of 
size n), with the following properties: 

1- PtPs = if t, s are different paths both of length n. 

2. zx = J2{Pt '■ t ends at A}. 

3. Pt® 1 = J2{Ps '■ s has length n + 1 and extends t} 

Let t be a path of length n and shape A and let /i be a Young diagram 
of size n + m. It follows that {pt (8) lm)Zfj. ^ precisely when there is a path 
on the Bratteli diagram from A to /i. It has been shown in that (in our 
notations) Tr(pt) = [Ai - A2 + 1], where [m] = (t"* - t"™)/(t - t"^) for any 
integer m, and where A is the endpoint of the path t. Observe that we get 
the same value for diagrams A and fi (of different sizes) that are in the same 
column in the Bratteli diagram. 

The idempotents pt were defined by recursive formulas in [W2], general- 
izing the formulas for the Jones- Wenzl idempotents in [Wl]. 

2.3 Specializations at non-roots of unity 

When r is not a proper root of unity, the Temperley Lieb algebras T„(r) are 
semi-simple complex algebras with the "same" structure as generic Temper- 
ley Lieb algebras. That is, T„(r) = 0;^TA(r), where Tx{t) is isomorphic to 
an fx-hj-f\ matrix algebra over C. The embeddings Tni^r) T„+i(r) are 
described by the Bratteli diagram as before. The idempotents pt, and the 
minimal central idempotents zx, in the generic algebras T„, are evaluable at 
r, and the evaluations Pt(r), resp. zx{t), satisfy analogous properties. 
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2.4 Specializations at roots of unity and evaluable idem- 
pot ents 

We require some terminology for discussing the case where r is a root of 
unity. Let g = r^, and suppose that g is a primitive i-th root of unity. We 
say that a Young diagram A is critical if w(A) := Ai — A2 + 1 is divisible by £. 
The m-th critical line on the Bratteli diagram for the generic Temperly Lieb 
algebra is the line containing the diagrams A with w{X) = ml. See Figure 
2.6. 

Say that two non-critical diagrams A and fi with the same number of 
boxes are reflections of one another in the m-th critical line if A 7^ /i and 
\w{X) — mi\ = \w{fi) — mi\ < i. (For example, with i = 3, [2,2] and [4] are 
reflections in the first critical line w{X) = 3.) 

For T a proper root of unity, the formulas for path idempotents in [Wl] 
and [W2] generally contain poles at r, i.e. the idempotents are not evaluable. 
However, suitable sums of path idempotents are evaluable. We will review 
some facts from [|GW|| about such evaluable sums . 

Suppose w{\) < i and t is a path of shape A which stays strictly to the 
left of the first critical line (in case w{X) < i), or hits the first critical line 
for the first time at A (in case w{X) = i); then pt is evaluable at r, and 
furthermore Tr(pj) = [w{X)]r = (r"'^ - t-'"^^'>)/{t - r"^). 

For each critical diagram A of size n, the minimal central idempotent z\ 
in T„ is evaluable at r. Furthermore, for each non-critical diagram A of size 
n, sua evaluable idempotent = '^pt ETn was defined in ||GW|| as follows: 
The summation goes over all paths t ending in A for which the last critical 
line hit by t is the one nearest to A to the left and over the paths obtained 
from such t by reflecting its part after the last critical line in the critical line 
(see Figure pT7| ). 

These idempotents have the following properties (which were shown in 
[GW]: 

1. {zx{t) : A critical } U {zj^^r) : fj, non-critical } is a partition of unity 
in T„(r); that is, the idempotents are mutually orthogonal and sum to 
the identity. 

2. zx{t) is a minimal central idempotent in T^It) if A is critical, and z^{t) 
is minimal central modulo the nilradical of T„ if A is not critical (see 

GW , Theorem 2.2 and Theorem 2.3). 
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Figure 2.6: Critical lines 



3. For A and non-critical, z^{r)Tn{r)zj-^{T) ^ only if A = /i, or if there 
is exactly one critical line between A and which reflects A to /x. If 
in this case u denotes the leftmost of the two diagrams A and then 
z^TnZ^ C Ti, (in the generic Temperley Lieb algebra). 

4. Let z!^^ = "^Pt, where the summation goes over all paths t which stay 
strictly to the left of the first critical line, and let 2;"*' = 1 — zl^^^. Then 
both zl^^^ and z"*' are evaluable; this is a direct consequence of the fact 
that zl^^ — z^, where the summation goes over diagrams A with n 
boxes with width w{X) < i. 

Proposition 2.1. The ideal of negligible morphisms in TL(t) is generated 
by the idempotent p\^^^{t) e r^_i(r). 
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Figure 2.7: A path and its reflected path. 



Proof. Let us first show that -z"*'(t) is in the ideal generated by P[^_i](t) for 
all n. This is clear for n < i, as = P[£-i] and 2;"*' = for n < £ — 1. 

Moreover, 2;"*' is a central idempotent in the maximum semisimple quo- 
tient of Tn, whose minimal central idempotents are the zj^ with w{\) > i. 
One checks pictorially that p[e_i]Z^ ^ for any such A (i.e. the path to [^ — 1] 
can be extended to a path t for which is a summand of z^. This proves 
our assertion in the maximum semisimple quotient of T„; it is well-known 
that in this case also the idempotent itself must be in the ideal generated by 
In particular, Hom(n, m)2;™^(T) -|- 2;™^(T)Hom(n, m) is also contained 
in this ideal. 

By [GW], Theorem 2.2 (c), for A a Young diagram of size n, with w{X) < 
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£, z^TnZ^{r) is a full matrix algebra, which moreover contains a minimal 
idempotent pt of trace Tr{pt) — [w{X)]t- ^ 0. Therefore 

z^r„z^(T)nNeg(n,n) = (0). 

Furthermore, z'!^^TnZ'!^^ {r) = zj^Tnzj^{T), by Fact 4 above, so 

C^7;C'(r)nNeg(n,n) = (0) 
as well. Now for x e Neg(n, n), one has z^^^ {t)xz^^^ {t) — 0, so 

a;eT„(r)C'(r)+zf(r)T„(r). 

We have shown that Neg(n, n) is contained in the ideal of TL(r) generated 
by for all n. That the same is true for Neg(m, n) with n ^ m follows 

from using the embeddings, and their left inverses, described at the end of 
Section 1.1. 

□ 



3 Ideals 

Proposition 3.1. Any proper ideal in TL (or in TL(r) ) is contained in the 
ideal of negligible morphisms. 

Proof. Let a G Hom(m, n). For all b G Hom(n,m), tr(6a) is in the intersec- 
tion of the ideal generated by a with the scaiars End(O). If a is not negligible, 
then the ideal generated by a contains an non-zero scalar, and therefore con- 
tains all morphisms. □ 

Corollciry 3.2. The categories TL andTL{T) forr not a proper root of unity 
have no non-zero proper ideals. 

Proof. There are no non-zero negligible morphisms in TL and in TL(t) for 
T not a proper root of unity. □ 

Theorem 3.3. Suppose that t is a proper root of unity. Then the negligible 
morphisms form the unique non-zero proper ideal in TL(t). 

Proof. Let J be a non-zero proper ideal in TL(r). By the embeddings dis- 
cussed at the end of Section 1.1, we can assume J fl r„ 7^ for some n. 

Now let a be a non-zero element of JnT„(T). Since {z\{t)} U {z^{t)} is 
a partition of unity in r„(T), one of the following conditions hold: 
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(a) b — az^{T) 7^ for some critical diagram /i. 

(b) b = zj^{T)azl^{T) 7^ for some non-critical diagram /x. 

(c) h = z^{r)azy{r) ^ for some pair A, A' of non-critical diagrams which 
are reflections of one another in a critical line. In this case, let /j, denote 
the leftmost of the two diagrams A, A'. 

In each of the three cases, one has b G e{r)Tn{r)f{T), where e, / are 
evaluable idempotents in T„ such that eT„/ C T^. Let a be a Young diagram 
on the first critical line of size n + m, such that there exists a path on the 
generic Bratteli diagram connecting /i and a. Then one has 

dime z^{T){e{T) ® l„)(T„(r) ® C lm){f{r) ® I J 

= dimc(t) )(T„ (g) C{t)lm){f <8) Im) 

= dimc(t) eTnf = dime e{T)Tn{T) f (t) 

where the first and last equalities result from the principle of constancy of 
dimension, and the second equality is because x i— Za{x ^ Im) is injective 
from Tn to T^. But then it follows that x i-^ Za{T){x (g) Im) is injective on 
e{T)Tn{T) f {t) . In particular (b ® lm)za is a non-zero element of J fl T^. 
Hence there exists c e Tq, such that g — c{b<^ im)za is an idempotent. After 
conjugating (and multiplying with P[e-i] ® Im, if necessary), we can assume 
to be a subidempotent of P[e-i] Cg) Im- But then Q-i+m/-i(fi') is a multiple 
of P[e^i], by Lemma 1.2, with the multiple equal to the rank of g in T^. This, 
together with Proposition 2.1, finishes the proof. □ 

It is easily seen that TL has a subcategory TL^*' whose objects consist of 
even numbers of points, and with the same morphisms between sets of even 
points as for TL. The evaluation TL^^(t) is defined in complete analogy to 
TL(r). 

Corollary 3.4. // is a proper root of unity of degree £ with i odd, the 
negligible morphisms form the unique non-zero proper ideal in TL^'". 

Proof. If £ is odd, P[e-i] is a morphism in TL*^^. The proof of the last theorem 
goes through word for word (one only needs to make sure that one stays 
within TL^^, which is easy to check). □ 
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